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Summary 


Experiments have been carried out in an attempt to detect linearly 
polarized radiation from the galactic spur at /11=30°. Five observational 
methods are discussed and the results of experiments using these methods 
are given. The frequencies used were 160, 178 and 408 Mc/s and the 
region of observation covered the more intense parts of the spur. All the 
experiments agree in fixing a low upper limit of a few per cent on linearly 
polarized radiation from the spur. 





1. Introduction.—It has been suggested by Tunmer (1) that the spur-shaped 
feature of radio emission extending from the galactic plane at /! = 30° and following 
almost a great circle across the sky is due to the spiral motion of relativistic 
electrons in the magnetic field along the local arm of the Galaxy. Such radiation 
is emitted with the electric vector in the instantaneous orbital plane of the 
electrons and Tunmer shows that if these electrons originate in the inner parts of 
the Galaxy and travel outwards along the arms, most of the radiation from the spur 
will be due to electrons with orbits which are nearly perpendicular to the magnetic 
field. Under these conditions, the radiation from the spur may be partially 
plane polarized and this paper is an account of some efforts to detect such polari- 
zation. 

Thomson (2), working at a frequency of 160 Mc/s using a 15° beam, has already 
reported an upper limit of 1 per cent for polarization in this region. The observa- 
tions described have the advantages resulting from the use of a smaller beam and a 
higher frequency. 

2. Experimental design considerations 

(i) Faraday rotation.—Several mechanisms associated with the Faraday 
effect increase the difficulty of observing polarization (3). For example, any 
disorder of the magnetic field gives components of the field along the line of sight, 
so that the plane of polarization is rotated by the Faraday effect. If the emitting 
regions in the line of sight are at different distances, there may be enough differential 
rotation of the vectors representing the contributions to the emission of the nearest 
and farthest regions to reduce considerably the observed polarization. It is to this 
mechanism that Tunmer attributes the low result of Thomson, because by taking 
values of o-1 cm~* for the electron density and 10~’ gauss for the longitudinal 
component of the field, she finds the rotation to be 0-03 radians per parsec at 
160 Mc/s. If the emitting region extends over 100-150 parsecs the rotation of the 
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plane of polarization across it may be 3-4 radians and the observed polarization will 
be small. 

Another mechanism which reduces the observed polarization is the differential 
rotation over the pass-band of the receiver. Both of these effects are reduced by 
observing at a higher frequency and the second by using a smaller bandwidth, 
but in both cases the sensitivity of the experiment is reduced. 

(ii) Spurious responses.—Aerial systems designed to detect plane polarized 
radiation may respond to unpolarized sources off the axis of symmetry as though 
these were polarized. It is necessary to take account of this effect, especially for 
paraboloidal reflectors of small focal ratio. 

(iii) Resolving power.—A narrow beam is desirable, first to avoid dilution 
of any polarized radiation from the spur (the width of which is some 5°) by 
unpolarized radiation from nearby regions and secondly to avoid the presence in 
the beam of regions with different planes of polarization. 

3. Observational methods. 

(i) Lonospheric rotation method.—The procedure here is to obtain drift curves 
of the total power received by an aerial with a small beamwidth over a period of 
several weeks, at a time of the year when the spur is in the beam at local dawn or 
dusk. The progressive displacement of sidereal time with respect to solar time 
causes the plane of polarization of any polarized radiation at the aerial to rotate 
as the electron density in the ionosphere changes. If a linearly polarized array 
is used, this rotation will change the total power received from day to day. 

At frequencies below about 200 Mc/s the method is useful since the change of 
rotation of the plane of polarization between day and night is large; at 160 Mc/s 
for instance, near the zenith, it is 7-8 radians. At higher frequencies the method 
is less satisfactory, the change of rotation at 400 Mc/s being only 1-1 $ radians. 

(ii) ‘Comb filter’ method.—To investigate the possibility of an appreciable 
rotation of the plane of polarization across the pass-band of the receiver a filter 
method can be used. 

In this method, the output of the aerial is split by a hybrid network into two 
parts which are fed into a phase-switching receiver through two cables. If the 
lengths of these are the same and the plane of the polarized signal does show 
several rotations across the pass-band, the output may be nearly zero. If the 
cables differ in length, the interference of the two parts of the signal produces a 
response which varies with frequency f as cos? (7d/c) f across the pass-band where d 
is the difference in the electrical lengths of the cables and c is the velocity of light. 
The phase-switching action of the receiver interchanges the maxima and minima 
of this pattern at the switching frequency. If the maxima of the frequency 
response and those produced by Faraday rotation of a polarized signal coincide, 
an output will be observed. 

The procedure is to observe with a given path difference, then with an additional 
quarter wavelength of path. This is repeated with other path differences, in 
steps of c/B where B is the bandwidth of the receiver. The useful range of path 
differences is limited at the lower end by the finite width and sharp edges of the 
I.F. response, but polarized radiation with a Faraday rotation across the pass- 
band smaller than this would be detected directly by method (i). 

(iii) The frequency-switching method.—This method is also designed to detect 
plane polarized radiation subject to appreciable Faraday rotation across the pass- 
band of a receiver as used in method (i). The aerial output is fed into a mixer, 
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in which the local oscillator frequency is switched at about 270¢/s by altering 
the potential across a ‘varactor’ diode in the resonant circuit, then through a 
narrow band filter into an I.F. amplifier, a switch frequency amplifier and a 
synchronous detector. If there is a Faraday rotation of any polarized com- 
ponent between the two narrow frequency bands accepted in this way, an output 
will be observed. A change in output will also be observed if there is a change 
of the spectral index of the radiation as different regions pass through the beam. 

(iv) Power difference method using crossed dipoles—The procedure is to 
measure the difference in the aerial temperature of crossed dipoles, for example 
at the focus of a paraboloidal reflector. A source of linearly polarized radiation 
will, in general, give an output at the receiver. If drift curves are obtained from 
the region under examination, with the dipoles in some arbitrary position and the 
observations are then repeated with the dipoles rotated through 45°, the 
magnitude and direction of the apparent polarization may then be calculated. 
An unpolarized source off the axis of the dish will also give an output because 
the polar diagrams of the dipoles in the FE and H planes are different. The 
spurious effects may be calculated by determining the polar diagram of the 
system for unpolarized radiation, using the Sun, and convolving this with a map 
of the brightness distribution over the sky. Subtraction then gives any polarized 
radiation that may be present. 

(v) Correlation method using crossed dipoles.—Correlation of the voltages 
induced in crossed dipoles at the focus of a paraboloid may be measured using a 
phase-switching receiver. As in (iv) an output is obtained both from linearly 
polarized sources and from unpolarized sources off the axis, although the spurious 
effects in the two methods are different and arise in different ways. Here, they 
are due to the cross-polarization terms of the polar diagrams of the dipoles (4). 
The procedure of observation and calculation of spurious effects is the same as 
that in method (iv). 

Alternatively, the ionosphere may be used to provide variations in the rotation 
of any polarized component of the radiation and the dipoles may be continuously 
rotated, so that all the information required is obtained during one run—this was 
the method used by Thomson. 

4. Observations and results 

In this section, the terms ‘ polarization temperature’ and ‘ degree of polari- 
zation’ are used. By the former is meant the difference in the aerial temperatures 
of two dipoles, one of which receives all the polarized radiation and the other none, 
and by the latter is meant this difference divided by the sum of the aerial 
temperatures. These definitions conform to normal usage. 

(i) Ionospheric rotation method.—Observations were carried out in 1958, 
March-April using a cylindrical parabolic aerial, at a frequency of 160 Mc/s. 
The beamwidth was 1°-2 in R.A. and 7°-7 in Dec., the bandwidth 2-5 Mc/s, and 
the region covered extended from R.A. 13" to 22" at Dec.+17°. The upper limit 
of 5 per cent on the degree of polarization was imposed by variations in receiver 
sensitivity. 

A similar experiment at 178 Mc/s was carried out in 1959 June, the bandwidth 
being 4 Mc/s, the beamwidth 13’-6 in R.A. and 4°-6 in Dec., and the region 
covered extending from R.A. 12" to 03" at Dec.+19°. An upper limit of 5 °K 
was put on the polarization temperature, corresponding to upper limits of the 
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degree of polarization between 0-6 per cent and 1°5 per cent over the different 
regions observed. 

(ii) The ‘comb filter’ method.—Observations were made in 1958 March at 
160 Mc/s with the same aerial constants as those given above and covered the 
region R.A. 10" to 21" at Dec.+17°. The cable lengths used corresponded to 
all possible Faraday rotations between 1 turn and 4 turns across the pass-band. 
Drifts of the zero level limited the accuracy of the experiment which gave an 
upper limit of about 5 per cent to linear polarization with the Faraday rotations 
indicated above. 

(iii) The frequency-switching method.—An experiment was carried out in 
1960 April, mainly as a test of the method, at a frequency of 178 Mc/s using the 
aerial of method (i). The switching took place between two bands of width 
250ke/s separated by 3 Mc/s. The region observed was at R.A. 00" to 24}, 
Dec.+19°. Drifts of zero level limited the accuracy, but an upper limit of 
5 per cent could again be placed on any linear polarization of the spur which rotates 
appreciably in the frequency interval given. Drifts due to variations in spectral 
index across the regions observed would be expected to be small compared with 
those observed. 

(iv) The power difference method.—Observations using a 7:5 m paraboloid at a 
frequency of 408 Mc/s were carried out during 1959 September. The beamwidth 
was 8°. The aerial was directed at the North Pole, then at the region R.A. 
17", Dec.+ 17°, which was allowed to drift through the beam, and finally at the 
North Pole again. Four runs were obtained with the same orientation of the 
dipoles; then the procedure was repeated for orientations at 45°, go° and 135° to 
the original. The spurious effects were calculated and gave a polarization 
temperature of about 0-5 °K. ‘The actual polarization temperature measured 
was 1 °K +1 °K, giving a degree of polarization in this region of less than 2 per cent. 

(v) The correlation method.—In 1959 October, the region of sky centred. at 
R.A. 17", Dec.+20° was followed with the same paraboloid as in the previous 
experiment for three hours after sunset. During this time, the orientation of the 
sky with respect to the dipoles changed by only 4°, so that spurious effects would 
remain nearly constant, whereas a rotation of about 60° was expected for any 
polarized component as the electron content of the ionosphere and the line of 
sight through it altered. An upper limit of 2 °K could be put on the polarization 
temperature, or about 2 per cent on the degree of polarization of the radiation 
from this region. 


TABLE I 
Per cent 
Method Frequency Bandwidth Beamwidth Region polarization 
a ) a °) less than 
(i) 159 Mc/s 2°5 Mc/s 12 7°%7 «61gh-22h, 0 4-19° 5 
(i) 178 4 13°6 4°6 125-032, +19° I 
(ii) 159 0°3-1'25 12 899° «= rob-arh, 44-17° 5 
(iii) 178 0°25 136 4°6 ob-2gh, +109° 5 
3 Mc/s apart 
(iv) 408 I 8° 17h +17° 2 
(v) 160 4 15° r7h-18h, +22° I 
(J. Thomson) 
(v) 408 I 8° 17h +20° 2 


(v) 404 4 8° r3%—170, +17° 402 
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In 1960 April-May drift curves were obtained of the region of sky 
R.A. 135-175, Dec. +17°, with the same paraboloid, at a frequency of 404 Mc/s, 
using two orientations of the dipoles. The spurious responses fixed the limit of 
detection of polarized radiation at 2 °K over the region as a whole (or a degree of 
polarization of 4 per cent near R.A. 13" to 2 per cent near R.A. 17"). No 
significant polarized radiation was detected. 

All these results are summarized in Table I, together with those of 
Thomson (2). 

5. Conclusions.—The experiments have shown that the degree of polarization 
in the bright region of the spur, at the lower frequencies of 160 and 178 Mc/s and 
at the higher frequency of 408 Mc/s, is very small and that depolarization across 
the bandwidth does not occur at the lower frequencies. 

It is possible to account for Thomson’s null result at 160 Mc/s on the basis of 
Tunmer’s mechanism if the usual values of the interstellar magnetic field and 
electron density are assumed, but it is difficult to explain the low results at 178 Mc/s, 
in view of the smaller beamwidth used, and those at 408 Mc/s, in view of the smaller 
interstellar Faraday rotation. Some alternative explanation of this galactic 
feature should therefore be sought. 

The theoretical implications of the observations will be discussed in more 
detail in a further paper. 


Mullard Radio Astronomy Observatory, 
Cavendish Laboratory, 
Cambridge: 
1960 November 4. 
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Summary 


It has been suggested that some apparently faint Fraunhofer lines may 
not be faint in the theoretical sense, but may be formed in low photospheric 
layers and ‘ veiled’ by a higher emitting atmosphere. 

This paper contains a theoretical investigation of this hypothesis. An 
exact solution is found for the residual intensity subject to certain standard 
though rather restrictive conditions. In the concluding section line profiles 
are drawn for veiled lines and for lines without veiling to show the effect of 
veiling on a line profile. 





1. Introduction.—In a letter to The Observatory C. de Jager and L. Neven (9) 
wrote: ‘‘It is well known to many workers in the field of stellar spectroscopy 
that some faint Fraunhofer lines need not necessarily be ‘faint’ in the theoretical 
sense, which means that they should not necessarily be situated on the 45° part 
of the curve of growth. There are lines that are caused by strong absorption in 
deep layers of a stellar atmosphere and which are veiled by the radiation from the 
non-selectively absorbing higher layers. Such lines may apparently be weak, 
but nevertheless behave as a strong or medium strong line, and show saturation 
effects since they are caused by strong absorption in the relevant parts of the 
stellar photosphere, where the line is principally formed.”’ 

In view of this suggestion, it seems desirable to use a model in which the star 
is surrounded by a continuously absorbing region, which will be called the atmo- 
sphere, and to compare the exact solution with the corresponding solution in the 
standard case in which there is no such atmosphere. ‘This comparison is effected 
by computing the profiles of veiled lines at various points of the disk in the two 
cases. More information about the centre to limb variation of the profiles of 
weak lines, thought to be veiled, will be required before adequate comparison 
with observation will be possible. 

The model used (see Fig. 1) is that of a plane parallel atmosphere of thickness 
x,, lying above a plane parallel photosphere of infinite extent. ‘Throughout both 
the atmosphere and the photosphere there is continuous emission according to a 
linear Planck function. Absorption lines are formed in the photosphere and the 
atmosphere is non-scattering. There is no radiation falling on the upper surface 
of the atmosphere. It is required to find the emergent intensity at any point of 
the line profile considered. 

In the first section of this paper an exact solution is found for the emergent 
intensity subject to certain standard conditions which will be enumerated later. 
The method used is one, developed by I. W. Busbridge, which was originally given 
by Ambartsumian for semi-infinite atmospheres. With this method it is 


* Received in original form 1960 August 3. 
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Fic. 1.—Diagram of the solar model used. The atmosphere is of thickness x, and the photosphere 
of infinite extent. In terms of optical depth the atmosphere has thickness 7, and in terms of the t 
variable it has thickness t,. I is used for the intensity in the atmosphere and I’ for that in the photo- 
sphere. 
comparatively simple to find an exact solution in terms of H functions, their 
moments, and integrals of expressions containing H functions. The coefficient « 
is introduced to allow for thermal emission and, in the numerical results, values 
for the emergent intensity are found in the case of lines formed by pure scattering 
(e=0) and pure absorption (€=1). 

A method of dealing with the variation of » (the ratio of the line absorption 
coefficient to the continuous absorption coefficient) is then developed in Section 
6 for the case of lines formed by pure absorption («= 1), and a formula capable 
of quite wide application derived*. 

Finally diagrams are drawn to show the effects of veiling on line profiles. 

2. Derivation of the Milne equation of the problem.—The problem is solved 
subject to certain conditions. ‘Throughout both the atmosphere and the photo- 
sphere it is assumed that 

(i) the variations of the continuous absorption coefficient x over the range 
of frequencies considered may be neglected ; 

(ii) the variation of the Planck function B(v, T) with the depth x below the 

surface of the atmosphere is such that it can be represented in the form 


B(v, T) =by + 5,7, 
where ad 
T= | Kp dx 


0 
and where the variations of the coefficients b) and 5, over the range of 
frequencies considered may be neglected. 


* Similar work has been done by Eddington (6) using a Milne-Eddington approximation and for 
a particular type of variation of 7 with optical depth. 
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Throughout the photosphere it is assumed, in addition, that 
(iii) the ratio », of the line scattering coefficient to « is independent of 7; 
(iv) the scattering in the line frequencies is isotropic, and also coherent, i.e. 
the frequency re-emitted within the line depends only on the radiation 
absorbed at the same frequency ; 
(v) the coefficient «, introduced to allow for thermal emission within the 
line, is independent of both v and r. 

Let r=7, when x=x, (see Fig. 1). In the atmosphere, (0o<7<7,), we 
denote the intensity in the line considered by J,(7, «), and in the photosphere, 
(7> 74), we denote it by J’,(7, ), a prime indicating the photosphere. In the 
atmosphere the intensity J,(7, 1), at an optical depth 7 and at an angle cosy 
with the outwards drawn normal, satisfies the transfer equation 


w(dldr)I,(r, 2) =1,(7, 2) — (by +647) (2.1) 
with the boundary conditions 

1,(0,—n)=0 (o<p<1), (22) 

L(t +h) =2 (ty +H) (o<u<1). (2.3) 


Now consider the radiation in the photosphere in the line frequencies. Here 
the intensity J’ (7, 1) satisfies the transfer equation 


mldlaa 0 Cr 6) = (140% 4) 4 —e)n [Pos ide! — (1401, (ba +047) 


(2.4) 
and the boundary condition 
I(t —4)=1' (ty — 1) (0<p<1). (2.5) 
If we write 
t,=(1 +7,)7, t=(1+7,)71, (2.6) 
A,= (1 +€n,)/(1 +75); n,=1/(1+7,), (2.7) 
and if we replace I’,(7, ~) by I’,(t,, u), then (2.4) becomes 
+1 
w(d/dt, (ty, 4) =1',(ty »)— (4 -A,) | V’,(tyy w!)du! —A, (bp + byn,f,) 
-1 
(2.8) 
and the boundary conditions (2.2), (2.3) and (2.5) become 
I,(0,—)=0 (0<p<1), (2.9) 
L(ty +h) =1'(ty +H) (o<p<1), — (2.10) 
L(ty—w)=1' (ty 1) (o<p<1). (2.11) 


We also need a condition limiting the rate of increase of J’ ,(t,, ») for large t,, We 
make this rate so small that the homogeneous Milne equation will have no non- 
null solution. (This restriction is in fact necessary to obtain a result which 
satisfies the physical conditions of the problem. For if we take the non-null 
solution we would obtain a source function which was exponentially large for 
large t,, and therefore a temperature which increased exponentially for large t,.) 

Using the condition (2.2) we first solve (2.1) for I,(7,—) and obtain the 
following expression 


1,(1,— 1) =by+5,(7—1) —(0p—byu)exp(—7/u) (0<m<t), (2.12) 











288 F. M. Hawkins Vol. 122 
and then the boundary condition (2.11) becomes 
T’ (ty, — 2) = bo + by (7, — 2) — (bo — yn) exp (— 74/1) 
= by +b,(n,t,— ph) — (bo— bn) exp(—2,t,/u) (0<p<1). 


(2.13) 
We now solve (2.8) subject to (2.13). Write 


+1 
B'(t.)=Hr-A)[ Flt Maa +24 (e+ amt) (2.14) 
-1 
then, dropping the subscript v for the time being, (2.8) is 
d 
SI(t,w)- =I, w= -=9'), (2.15) 
be Ke 
giving, for o<p<1, 
I'(t,+m)exp (—thu) = ["8'(=)exp (—s/u)n-* de. (2.16) 


(This assumes I'(t, +) exp(—t/u)—->0 as t 00). 
When p <0, write —y for uw in (2.15). Then, foro<p<1, 


I'(t—w)exp (tH) ~1 (ty —n)exp (hls) [ 8G) exp (x mud, 


(2.17) 
where /'(t,,—) is given by (2.13). From (2.16), (2.17) and (2.14) we obtain 


B'()=Ha—2) [TC w Yd + Abo bat) 
=r -a{ | 1 +p! )dp! + | 1 " Hn’) +A(by + 5,nt) 
=4(0—2)[ aw [0 G)exp [- (e—H/u' ur ae 
+402) ae’ 8'(x)exp [— (t=) Jw) 
+ H1—2)[ (ty —H') exp [— (tf )/ pd’ + A(by + byt 
ie. O(t) =4(1 a | “o'(e)de | “exp [—u(x—2)] udu 


+4(1- af 3’ (x) de “exp [—u(t—x)]u-! du+ B(t), (2.18) 
where u=1/p’ and 
Be) = 41-9) [ E(t —H exp [= (t= f)/o' de! +A(by + Dyn). (2.19) 
Equation (2.18) is 


8’ (t)= (1 —A) | “8'(2) Ex(e t)dx+3(1—A) | j 8’ (x) E,(t—x)dx + B(t), 
; (2.20) 
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14) 
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where 
E,(s)= | “x” exp (—sx) dx. 
1 
The right hand side of (2.20) is not quite 


(x—A) A{S'(x)}+ Bie), 


where 
A, tI} =4 | “IE, (t= at 


Write t=t,+y; x=t,+4u, so that u=y when x=2. 
Then (2.20) is 


8'(t,+y)=H(1—-A) | “8 (+) By(u—y)du 


+4(1—2) | "8 (ty +u) E,(y—u)du + B(t, +y) 
0 
and this is the Milne equation 
B'(t4,+y)=(1-A) A, {TF (t, +u)} + Bit, +y), (2.21) 
where, from (2.19) and (2.13), 


B(t)=3(1- a) fe +b, (nt, — pw’) — (by — by") exp (—nt,/u')}exp [— (t—t))/m’]dy’ 
+A(bo + 5, nt) 
= 3(1—A{(bo + bynt,) Ey(t — t,) — 6,5 (t — t,) — bo E g(t — t, + nt) 
+ b,E,(t—t,+nt,)}+A(bp + ,nt). (2.22) 
Define B,(y) by 
By) =B(, +9). (2.23) 
Then on putting t=t,+y in (2.22) we obtain 
By(y) = 4(1—A){ (do + bynt, )/Eo(y) — by E3(y) — boa (y + nty) + 5,Es(y + nty)} 
+A(by + bynt, + b,ny). (2.24) 
On writing B’(t, +y)=F(y), (2.25) 


(2.21) becomes 
F(y)=(1—-A)A,{F(u)} + By(y). (2.26) 


From (2.16) 

(ty +n)= [ “8'(e)exp [-@—4)/u] uta, 
and on writing ¢, +y=x this becomes 

(ty +H)= [8 +y)exp [-ylale ray, (2.27) 
i.e. 


I(ty +n) = | -F(y)exp [—ylu]u-*dy. (2.28) 


The next section will be devoted to finding results which are needed to find 
I'(t,, +) when F() is a solution of (2.26). 
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3. Preliminary mathematics.—We begin by solving the equation 
(1-wA),{F,(t)}= E, (7 +4) (3-1) 
in which w=1—A, a>o and n=2, 3... . First we consider the auxiliary 
equation [see (2), section 3] 
(1 -wA), (F(t, 8)}=exp (—87). (3.2) 


We can multiply by exp (— a5)d65/5” and integrate with respect to 6 over (1, 00) and 
this gives 


(1 -wA), | “exp (— a8) F(t,8)3-" as} 
1 
=| “exp [—8(7-+a)]8-"d3=E,(7-+a). 
1 
Thus 
F(t)= | “exp (— a8) F(t,8)8-" d8 (3.3) 
1 

is a solution of (3.1). It is known [see (2), theorem iv] that 

Gy, {F(t,8)}= | “F(ty8)exp (—t[u) ude 


ee 


eh (3-4) 


where H(,) is the solution of 
1 
H(y)=1+ Hue ae dé (3.5) 


which is continuous for o< <1 [see (4), Chapter 2]. Hence, by (3.3) 


Ii(ty +n) = | “Fy(y)exp (=|) w-2dy 
= | “exp (—ylbe) ut dy | “F(y,8)exp (~a8)5-* db 


= | “exp (—a8)3-" ab [ “#(,8)exp (—y/u)u-*dy 


and therefore, by (3.4), 


Ts (ty +n) = j Me exp (—a8)8-n ds. 
On putting = 1/£, this becomes 
Tally +H) =H) | — H(E)exp (—a/é) dé. (3.6) 


In (2), theorem v (with y= 1 —A), it is shown that the Laplace transform of the 


solution of 
F,(y)=(1-A)Ay {F2(t)} + a7 + (3-7) 


. 


T(ty +4) = | ~Fa(y)exp (—y/h) 0-4 dy 


= H(w AM? (dope + ay) + bad (1 —A)ay}, (3-8) 
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where 
1 
tty = | H(é) & dé. (3.9) 
0 


4. The solution for the emergent intensity.—The results of Section 3 can now be 
used to find the Laplace transform /'(t,, +) of the solution of (2.26) with B,(y) 
given by (2.24). As the operator A, is linear, this solution is 


I (ty +H) =A Oy thant) HH) [ pene 
—b,H(u) er F— H(é)ae 
— by H(u) ie: “pa Hle)exp (—nty/é) a 
+b, H(u){ HG) exp (—neylé) ab} 
+ H(p) {A¥? (b, np + by + bynt,) + $(1 —A) ab,n}. (4.1) 
By using (3.5) and the relations 
m=1—A, 3$(1—-A)a=1—Al? 


[see (2), equation (4.15)] the following expressions are found which simplify 
equation (4.1): 


H(u) fiz zi H() a= H(u) es ete= —H HE) a 


= H(t) a%— =a) [H(u)—1] 
=2(1—A)> [1-2 A(u)], 


H(u) ie A H(6) d= Hu) f. SEAS HCG) at 
ar apeere [1A H(y)]. 


On substituting these in (4.1) and writing 7,=mt, and 


Stray) = 401 -AH(u) {6 =£— H(@) exp (— rE) df 





—b,[ — H(eyexp(—n/é)at}, (4.2) 
we obtain 
I'(7y) + 2) = (bo + 4474) [1 — V7 (u)] — 4(1 —A)by, (4) 
+ byp [1 —A¥? H(w)] + A¥? (bymp + bo + 6,7) H(u) 
+ (1 —A)ob,nH(u) — S(74, 4); 


i.e. 


V(r) +1) = (by + by7y +52) — XY? yy (1 —n) H(u) 
— (1A) byey (1 —n) H(u) — Sry H)- (4-3) 
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The solution of (2.1) subject to (2.3) is 


I(o, +) =exp (— 74/4) 1'(7, +H) +o [exp (—7/)(b) +5,7) dr 


=exp (—7,/H)1'(7, +4) +59 +p 
— exp (— 74/-)(b9 + 5,7, + by) (4-4) 


where I'(7,, +) is given by (4.3). 

For emphasis J(0, +) will be written as J, |, (0, +), and the corresponding 
quantity in the standard case in which there is no non-selectively absorbing 
atmosphere will be written as J, ,(0, +). Theintensity in the continuum (which 
is the same in both cases) will be denoted by I*(0, +p). Let r, , and rp , be 
respectively the percentage residual intensities for the veiled line and for the line 
in the standard case. 


Then pap a 10M» (0, +H) —Fo( +H)] (4.5) 
Ty 0,v I*(o, +p) . 


The emergent intensity in the standard case in which there is no non-selectively 
absorbing atmosphere is given by* 


To,, (0, +) = H(u) {A¥? (by + dyn) + 3(1 —A) ob ,n}. 
Hence, from (4.5), (4.4) and (4.3), 


ae te) fp [7,,,»—%o, »] =exp (— 7/1 ){Do + by 7, + bys —AY7b, p(t — 2) H(y) 


10 
— $(1—A)byaq (1 — 2) H() — S(7y, w)} +49 + Op 
—exp (— 74/)(bo + 417, + by) 
—8(by + bynp.)H(u) ~ 4(1 —A)aybyn A(q) 
= —exp (—7/n){A"? d,n(1 —2) (yp) 
+ §(1—A)b,a4(1 —n)H(u)} +p + by 
— 2 (by + bynp) H(u) 
— $(1 —A)a bn (1) — exp (— 74/H) S(7y 1). 
On substituting for S(7,,) from (4.2) we get 
P+ H) Ey, ,— roc] = —exp (— rule) H(H) (bo + bye) + $04(1—A)a} 


100 
— H(p)[1 — exp (— 73/1) A"? (do + bymp) 
+ 4(1—A) dyna} + by +b yp 


~ }(1-A) exp (— nile (1) bo | oH) exp (— 1,/é)dé 





~b, {5 H(@)exp (— nie ae} (4-6) 
where J*(0, +)=b) +b, [see (5)]T- 


* See (5) Section 84, equation (66). This is also obtained from (4. 1) on putting t,;=o. 
+ The value of I*(o, +) can also be obtained by letting 7,->00 in (4.4), viz. I*(0, +4)=bo +b 
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In the case of lines formed by pure absorption, (« = 1), (4.6) can be reduced toa 
much simpler expression. 

On putting «=1 in (2.7) this gives \=1 and, from (3.5) witha =1—A=o, we 
see that 


H(y)= t. 
Hence (4.6) reduces to 
I ote) [r,,,»—%o, v] = 4m [1 — exp (—74/H)] [1-2], (4.7) 


and the expression for the percentage residual intensity in the veiled line becomes 


1*(,+n), 
100 


1, v= Oo t+ Opt — (1—n) exp (— 7,/)}. (4.8) 


5. Numerical results.—A plane parallel model is not a good representation 
near the limb and numerical results are therefore obtained only for o-2<p<1. 
However it is interesting to note that in every case 


I*(0, +“) 
oO 


— > (b9 + byp),, 20 as | i +0. 


Ts, 9 
Since this is ]*(0, + 4),,_ 9, the lines fade out completely at the limb. 

As any line is formed by a mixture of absorption and scattering it was decided 
to consider two extreme cases, one in which the lines are formed by pure scattering 
(e=o), and the other in which the lines are formed by pure absorption (e=1). In 
both these cases values of 7, ,,—7,, and r,, have been calculated for various 
values of 7, and for two thicknesses of the atmosphere, using the new tables of H 
functions prepared by Stibbs and Weir [see (10)]. In the former case the elec- 
tronic computer at the University Computing Laboratory was used in obtaining 
numerical results and my thanks are due to the Director for allowing use of these 
facilities. The results are given in Tables I-IV. 


Tasle I 


The veiling effect of the atmosphere (€=0) 


=” I o°9 o8 o"7 06 o'5 "4 03 o-2 
v/] \ Torn, »— 70, 
I 6°11 6°54 7°03 7°61 8-31 9°16 10°25 11°68 13°68 


4 10°89) )=—ssarI'79—s—«*a286— ss 14'16—sd15*BO=—«a17"QI_-~Ss 20°74 24°74. Ss 30°80 

9 13°20 14°37 15°79 17°54 19°74 22°62 26°53 32°14 40°70 
19 14°86 16°24 17°92 20°00 22°64 26°09 30°81 37°58 47°98 
99 17°20 18°88 20°92 23°45 2666 30°88 36°64 44:94 57°66 
399 18:23 20°02 22°20 24°90 2834 32°84 39°00 47°87 61°46 


Tors, v—To v 
I 8-65 9°20 9°82 10°53 11°37 12°94 «=: 13°§t 14°90 16°49 
4 15°37. 16°54 17°92 19°57 21°59 2411 27°33 31°59 37°24 
9 18-60 20°15 21°98 24°21 26:96 30°45 34°98 41°06 49°28 
19 20°92 22°74 24°03. 27°59, 30°90 F351 = g0'62 48:04 «= 58-173 
99 24°18 26°40 29°07 32°32 36:37 41°55 4833 57°48 69°93 
399 25°61 27°98 30°84 34°32 38°66 44:19 51°46 61:25 74°56 
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73°71 
53°70 
44°24 
37°50 
27°99 
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76°25 
58-18 
49°64 
43°56 
34°97 
3119 


5°44 
8-70 
9°79 
10°33 
10°79 
10°85 


7°78 
12°44 
14°00 
14°77 
15°40 
15°51 


75°44 
60°70 
55°79 
53°33 
51°39 
51°00 


77°78 
64°44 
60-00 
57°77 
56-00 
55°66 
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TaBLe II 
The centre to limb variation of the veiled lines (€=0) 


o'9 08 o"7 06 o'5 o*4 
To-2,¥ 
7510 7666 7843 8047 82°83 85-61 
55°59 57°74 60:20 63:10 66°57 70°86 
46°14 4835 50°93 53°90 57°80 62°65 
39°36 = 41°55) 44°15) 47°33. 51°32 = 56°58 
29°81 863198 §=— 34°65 37°99 42°33 48°19 
25°64 2736 30°59 34°06 38:58 44°76 
To°3, v 
7776 79°45 81°35 «= 83°53. 86°01 = 88-87 
60°34 62°80 65°61 68-89 72°77 77°45 
51°92 54°54 57°60 61°22 65°63 7110 
45°86 48°56 51°74 55°59 © 60°34 66-39 
37°33. 40°13. 43°52, 47°70 53°00 59°88 
33°60 36:50 gororn «44°38 = 49°93, 5:7°22 
Tase III 
The veiling effect of the atmosphere (€=1) 
o'9 08 o'7 06 05 o"4 
Tora, vo, v 

5°72 6:03 6°37 6°71 7°07 7°38 

9°16 9°66 10°19 10°74 11°31 11°80 
10°30 10°86 11°46 12°09 12°71 13°28 
10°87 11°46 12°09 12°76 13°42 1401 
11°33 1I'94 12°60 13°29 13°99 14°61 
11°42 12°04 12°70 13°40 14°09 14°72 

To-3,v—To, v 

8-14 8°53 8-93 9°32 9°67 9°89 
13°04 135 1429 1491 15°48 15°83 
14°66 15°33 16°05 16°77 17°40 17°81 
15°47 16-20 16°96 17°71 18°37 18°79 
16°12 16°89 17°68 18°45 19°14 19°59 
16°25 17°02 17°81 18°59 19°29 19°74 

TaBLeE IV 

The centre to limb variation of the veiled lines (€=1) 

o'9) 08 o'7 06 0°5 o"4 

To-2, 
77°00 =- 78°76 —S 80°76 = 83°03 85°64 «= «8888-63 
63:20 66:02 69:21 72°85 77°02 81°80 
58:60 8 61°77, 65°36 = 69°46 7414 79°53 
56°30 59°64 63:43 67°76 72°71 78:39 
54°46 57°94 61°89 66:40 71°56 77°48 
5412 5763 6161 6615 71°34 77°31 
Tors, v 

79°42 81:26 83:32 §©=685°64 «Ss (8824s “14 
67:08 70°01 73°31 77°02 81-19 85°83 
62°96 66:24 69°95 74:14 78°83 84:06 
60°90 }«=-_s«d64"38 Ss“ 68-30—s-772°71 77°66 83°17 
59°25 62:89 66:97 71°56 76°71 82-46 
58:95 6261 66°72 71°34 76°54 82°33 


o*3 


88-94 
76°38 
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56°56 
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92°16 
83°23 
78-12 
74°31 
69°10 
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15°07 


9°81 
15°69 
17°65 
18-63 
19°42 
19°57 
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14°34 
16°14 
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95°76 
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91°93 
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6. Variation of n, with optical depth.—Since the objection may be raised that 
the variation of », with optical depth has been neglected, a method will now be 
developed which will take into account this variation in the case of lines formed 
by pure absorption («=1). 

The same notation will be used as in Section 2, except that 7, will be denoted 
by 7,(7) as it is a function of optical depth. The problem will be solved subject 
to conditions (i), (ii), (iv) and (v) of Section 2. 

In this case equation (2.4) becomes (with «= 1) 


w(d/dr)l’, (7, n)= [1 +9,(7) (7, #) — [1 +0,(7)] [bo + 6,7] (6.1) 
and from (2.1) and (2.10) the expression for the emergent intensity is 
T,,, (0, +h) =1'(7, +H) exp (— 71/H) +49 + bm 
— (bo +6,7, +n) exp(—74/u) (0<p<1). (6.2) 
Equation (6.1) has the integrating factor 


exp{— f(r +n,(0]uath, (6.3) 
and on writing 
n(r=afrtn(yh x)= inde (6.4) 
(6.3) becomes 
exp [-9(7)/a]} (6.5) 


On multiplying (6.1) by exp[—y(r7)/u] and integrating with respect to 7 over 
(7, 00), the following expression is obtained when o<p<1: 


al'(rs + wdexp[—y(7)]= [bo +bat) Lx +n, (0) exp [—9(e)/u] at 


(6.6) 
[This assumes J’(7, +) does not increase exponentially for large r.] 
Since 7,(r)>0, therefore 
1/n,(r)=1+,(7)>1 
and hence, for r> 7}, 
wtr)= [fn (t)) de> 71 (6.7) 
Thus y(7,)=0 and y(r)—> 0 as 7 00. On putting r=7, in (6.6) this becomes 
pl'(r, +n)=| (bo + byt) [1 +7,(¢)] exp [—y(4)/H] at. (6.8) 
This expression can be simplified since, from (6.4), 
(d/dr) y(r)=1/n,(r)=1+7,(7), (6.9) 


and therefore 
(d/dr) exp {—y(r)/u}= —w [1 + 9,(7)] exp [—y(7)/+]. (6.10) 
Also, from (6.7), 


(bo +67) exp[—y(7)/m]>0 (7 > &). 
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Hence, on integrating (6.8) by parts, we obtain 
pl'(7, +u)= [ — (dy + b,t) exp {—y(t)/p}] 


+ pb, | “exp [—p(0)/u] dt 


= ps (by + by 74) + wy | “exp [~y(0/q] at (6.11) 
Then from (6.2) the emergent intensity is given by 
Ty (0, +) = (9 + 5; 71) exp (— 7,/H) 
+byexp(—n/u)| “exp [—y(e)/u] at 
+ by + bp — (by) + by 7, + by) exp (— 74/1) 
= by + bp [1 — exp (— 74/4)] 


+b, exp (—74/u) | “exp [—y(¢)/1] dt. (6.12) 


The expression (6.12) is an exact solution for the emergent intensity in lines formed 
by pure absorption, and no particular law of variation of 7, with optical depth is 
presupposed. 

Numerical values for the percentage residual intensity were obtained in the 
case in which the photosphere is divided into two layers, in the upper of which 
n,(7) was a non-zero constant and in the lower 7,(7) was zero. The values chosen 
for 7,(7) in the upper layer were some of those used in the first part of this paper, 
viz. 1, 4, 9 and 19, and for these the optical thicknesses of the layer were taken to be 
0°5, 0°2, o'r and 0°05 respectively. The values 0-2 and 0-3 were again used for 7,. 
By comparing the results thus obtained with similar ones obtained by I. W. 
Busbridge for 7, =0 (see (3)), values for the veiling effect of the atmosphere were 
also obtained. Since, however, all my results for this case are similar to those 
obtained previously (see Section 5), it seems unnecessary to insert tables of them. 

7. Conclusions.—In the previous sections formulae have been obtained which 
can be used to find the percentage residual intensity in certain cases. It now 
remains to compare the results, so obtained, with relevant observational data and 
to show how a veiling atmosphere would affect a line profile*. In Fig. 2 the varia- 
tion of the percentage residual intensity for 7=1 is compared with observational 
data for weak lines obtained by M. G. Adam for faint Fraunhofer lines (1). The 
full line shows the variation of the percentage residual intensity from centre to 
limb for «=o and 7, = 0°3 obtained from the formulae of Section 4. The scattered 
points show the values for the central intensity (expressed as a percentage) 
obtained from Miss Adam’s work. ‘The filled circles are for the line A 5102-975, 
the triangles for 45145-470 and the crosses for 1 5147-484 (these lines have res- 
pectively Rowland Intensities of 1, 0, 0), where for each line the mean of the results 
from three plates is shown. ‘These lines are weak and are chosen because they 
give points quite close to the full curve at »= 1. 


* For observational work on line profiles see, for example, (7) and (8). Computations have been 
made by de Jager and Neven (ga) for solar lines of N1 and O1. Working from curves of growth 
for certain lines of N1 and O1, computed with the aid of the Pecker saturation function, they showed 
that large errors are introduced by the usual weak line approximation. 
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In this paper the parameter 7, has been used to denote the ratio of the line 
absorption coefficient to the continuous absorption coefficient. For a line 
produced by a certain atom, therefore, 7, is the ratio of the atomic absorption 


Intensity 


100 T T T T T 














70 ! | 4 l | i | 
] 0-8 0-6 y 0-4 o2 





Fic. 2.—Comparison of the variation from centre to limb of the residual intensity, obtained 
theoretically, with M. G. Adam’s observational results. The ordinates represent residual intensity 
(expressed as a percentage). 


coefficient to the continuous absorption coefficient. As we can take the variation 
with frequency of the continuous absorption coefficient over a particular line to be 
negligible, the coefficient », varies as the atomic absorption coefficient over the 
line. In the centre of the line (Doppler core), therefore, 7, is proportional to 
exp (— AA?/A)?,,) and in the wings (damping wing) 7, is proportional to 1/A)?, 
where Ad is the distance from the centre of the line in units of wavelength and 
AX, is the Doppler half width (i.e. AA; is the distance from the centre of the line 
to the point at which the intensity has e~! of its value at the centre of the line). 

In Fig. 3 the line profiles were drawn using 7 = 399 exp (—v”) for the central 
part of the line and 7 = g/v? for the wings, and those in Fig. 4 using 7 = 19 exp (—v?) 
for the central part of the line and »=2-25/v" for the wings, where v=Ad/AA>. 
In each case 7 is taken to be independent of depth in the photosphere. These 
figures show how the line profile of a typical Fraunhofer line varies from centre to 
limb under certain circumstances. In each figure, the line profile at the centre of 
the Sun (u=1) and at the limb (~=0-2) are shown. Fig. 3 shows the case in 
which the lines are formed by pure scattering (e=0). In Fig. 3(a) there is a 


22* 
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Fic. 3.—Line profiles at the centre (u=1) and at the limb (u=o0-2) of a typical Fraunhofer line 
for the case in which the line is formed by pure scattering (e=0). In Fig. 3(a) there is a veiling atmo- 
sphere of optical thickness 0-3 and in Fig. 3(b) there is no such atmosphere. In each case the distance 
from the centre of the line is measured in units of 4{4Ap, and the ordinates represent percentage 
residual intensity. 


veiling atmosphere of optical thickness 0-3 and in Fig. 3 (b), which is shown for 
comparison, there is no such veiling atmosphere. Fig. 4 shows the case in which 
the lines are formed by pure absorption (e=1). In Fig. 4 (a) there is a veiling 
atmosphere of optical thickness 0-3 and in Fig. 4 (b), which is shown for comparison, 
there is no such veiling atmosphere. 
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percentage residual intensity. 
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These figures show the behaviour of the line profiles in two extreme cases, viz. 
lines formed by pure absorption and lines formed by pure scattering. As there is 
some similarity between Figs. 3 (a) and 4 (a), it would not be difficult to visualize 
the variation of the line profile of a veiled line which is formed by a mixture of 
scattering and absorption. 

In conclusion, I should like to thank Dr I. W. Busbridge for all her help and 
encouragement, and also Professor H. H. Plaskett who made several very helpful 
suggestions, Professor D. W. N. Stibbs who advised me about the best method 
of integrating expressions containing H functions, Dr W. R. Hindmarsh for very 
useful discussions on the physics of absorption lines, and the Department of 
Scientific and Industrial Research for awards, during the tenure of which this 
work has been carried out. 

University Observatory, 

Oxford: 

1960 September. 


References 


(1) M. G. Adam, M.N., 98, 112, 1937. 

(z) I. W. Busbridge, Quart. 7. (Oxford), (2), 6, 218, 1955. 

(3) I. W. Busbridge, M.N., 116, 304, 1956. 

(4) I. W. Busbridge, Mathematics of Radiative Transfer (Cambridge, 1960). 

(5) S. Chandrasekhar, Radiative Transfer (Oxford, 1950). 

(6) A. S. Eddington, M.N., 89, 620, 1929. 

(7) G. E. Hale and W. S. Adams, Ap. 7., 25, 300, 1907. 

(8) J. Houtgast, The Variations in the Profiles of Strong Fraunhofer Lines along a Radius of 
the Solar Disc, thesis, Utrecht, 1942. 

(9) C. de Jager and L. Neven, Observatory, 79, 102, 1959. 

(ga) L. Neven, Comm. Obs. Roy. de Belgique, 157, 112, 1959. 

(10) D. W. N. Stibbs and R. E. Weir, M.N., 119, 512, 1959. 








es re 


—-— -—_—~ 


oo == -& fF &§ 4 Ht © ff 





122 


VIZ. 
re is 
alize 
e of 


and 
pful 
chod 
very 
t of 
this 


us of 





A STUDY OF SUNSPOT VELOCITY FIELDS USING A 
MAGNETICALLY UNDISTURBED LINE 


Fj. Holmes 


(Communicated by the Director, University Observatory, Oxford) 
(Received 1960 November 28) 


Summary 


The relation between magnetic lines of force and material motion is a 
matter of considerable importance in sunspot theory, and this investigation 
of sunspot velocity fields uses a line for which the Landé splitting factor 
g=o. (Fer A 5576-101). Doppler displacements were measured at over 
100 points in the region of both a large and small sunspot. At the high 
dispersion value of 0-17 A/mm, measurement difficulties are presented 
by the broad solar lines, but a check on visual accuracy was provided by 
a new photo-electric method of measurement. The measured velocities 
were corrected for observer’s motion and solar rotation and reduced 
to velocity components, with origin at the spot centre. Using Kinman’s 
conclusion that the spot flow is wholly radial, we have obtained fair agreement 
for the maximum radial velocity in the penumbra of the small spot but found a 
value much smaller than that predicted by Kinman for the large spot. 

The differences from the earlier measures have been attributed to the fact 
that we are now able to measure the true Doppler shifts, freed from the 
complicating effects of Zeeman splitting, and also to the “ line-flare ’’, which 
appears in the penumbral region of high dispersion spectra, but was not 
resolved in the earlier work on Evershed effect. 





Introduction.—Motion in the neighbourhood of sunspots was first investigated 
by Evershed in 1901, and its general form is now well established (1, 5, 6, 7, 8, 
g, 18, 19). Detailed measurements were made by Kinman at Oxford in 1952 
(10, 11) and again by Servajean at Meudon in 1958 (17). What has been neglec- 
ted to date, however, is the distinction between line shifts due to velocity along the 
line of sight and shifts which have their origin in the line splitting produced by the 
magnetic field of the sunspot. ‘The sunspot velocity fields and magnetic fields 
are both observed to change in a more or less regular fashion across the penumbra, 
and the relation between magnetic lines of force and material motion is a matter of 
considerable importance in sunspot theory. Plate 6 (a) shows the splitting of the 
Fe line 46302 across a sunspot. This clearly shows that the partial suppression 
of one or other of thetwo o components of the pattern due to instrumental polariza- 
tion can easily lead to an apparent line displacement. A line such as A 6302 is an 
extreme case and would not be used for velocity measures, but such effects exist 
on a smaller scale for all lines, except those for which the Zeeman splitting factor 
iszero. Itis only by using such a magnetically unaffected line that we can be sure 
that the velocities we measure are free from magnetic field effects*. 

For this reason we have selected the Fe line 4 5576-101 for our present study 
of sunspot velocity fields. This is a line first suggested by Von Kliiber (20), for 

*Servajean’s work (Ann. d’ Astrophys., 24, 1, 1961), using the line \ 5691-508 (Fe Ni) 
for which g=o, has been published while the present paper was in the press. 
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which the Landé splitting factor g=o. But we have now to contend with the 
observational difficulties presented by such a choice. There are no terrestrial 
standards in this region, such as the oxygen lines which appear in Plate 6 (a), or the 
water vapour lines used by Kinman, but the difficulty has been overcome here by 
introducing absorption lines due to iodine vapour into the solar spectrum. These 
could not be introduced directly into the sunspot spectrum because the iodine line 
density is too high but, instead, two exposures were made in rapid succession; 
the first for the spot spectrum alone and the second for the spot spectrum plus 
iodine. The iodine lines could then be used to show slit curvature and possible 
slit inclination, and consequently these effects could be allowed for in the solar 
exposure. We then assume that at the top and bottom of the slit, far removed 
from the spot, the line shows only the normal photospheric velocity shift due to 
solar rotation, plus a possible limb effect, if we take the effect of the photospheric 
velocity fields to be uniform across the slit. After calculating and allowing for 
such effects at selected points across the slit, we have the true Evershed effect. 

Observations have so far been made on one small and one large spot. Central 
meridian passages (kindly supplied by the Royal Greenwich Observatory) were 
September 11.30 and September 3.25 respectively. The observational material, 
together with its measurement and reduction, are described in the first two sections 
of the paper and we conclude with a short discussion of the velocity fields now 
obtained. 

1. Observations and measurements.—The sunspot spectra for the investigation 
were obtained by Professor Plaskett in September 1959. The image radius py as 
given by the new solar telescope was 165mm and the slit height of 25 mm 
corresponds to approximately 145 seconds of arc or approximately 105 500 km on 
the solar surface at the disk centre. The Babcock grating was used in the 5th 
order with the 40 ft spectrograph, resulting in a linear dispersion of 0-17 A/mm. 
In each exposure, the slit passed through the spot centre and details of the mea- 
sured plates are given in Table I. In this table, p=radial distance of the spo 
from the disk centre and B, L are the heliographic latitude and longitude of the 
spots. 


TABLE I 
Measured plates 
Measured line A 5576-101 z°F°—e®D Fel. 
Slit height =25 mm. 
Slit width =0-045 mm. 
Dispersion = 0°17 A/mm. 


Date Exposure Umbral 
iota 1959 U-T. time . . P/PO diameter 
E22 16202M 368 40sec iis = E 46” 
E23 Sept. 7 165 72M 358 60 sec ¥57 1S 54 0°793 34000 km. 
E28 08) 41M 578 30 sec , * 
E29 oe.) A 49™388 31sec — ” 
—18° 10’ 274° 50’ 12000 km. 
E36 08) 28M 168 30 sec ’ 
E37 Sept. ro og) 17M 068 30 sec — 


Measurements were made using the Hilger micrometer with a transverse 
motion for the plate stage (16). The field of view, perpendicular to the dispersion, 
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(a) Zeeman splitting in Fe A 6302. 
(b) “ Line-flare ”’ in large and small spots. 


J]. Holmes, A study of sunspot velocity fields using a magnetically undisturbed line. 
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was restricted to 0-3 mm by a diaphragm in the microscope eyepiece, so that this 
width of spectrum can be measured at any desired height. Since the sunspot 
spectra contain no atmospheric lines or lines of zero velocity, a ruled glass plate 
was used to provide a fiduciary line. This was permanently attached to the 
microscope stage and check settings were made on the ruled line throughout all 
the micrometer measurements of sunspot and iodine lines. Provided the sunspot 
spectra and iodine spectra are closely aligned at the same angle with respect to the 
fiduciary line, the iodine line readings give a correction curve as a function of slit 
height, by means of which all instrumental shifts may be removed from the 
observed displacements of the solar line. The six iodine correction curves were 
identical within the errors of measurement and were combined to give a single 
curve by means of a least squares solution. 

All spectra were measured in the direct and reversed directions and mean 
readings obtained. ‘The sunspot spectra were measured at over 100 points over 
the 25 mm slit height at intervals as small as 0-1 mm in the penumbral region, and 
wider intervals in the outer parts, where we suspect the Evershed velocities are 
changing less rapidly. A complete measurement of the plates was made in 
February and the whole set of readings repeated in April as a check, final readings 
being taken as the mean of these two sets. Although the high dispersion is a 
great advantage in the work, the solar lines are now so broad as to make visual 
measurements a difficult matter, especially over the restricted heightourdiaphragm 
imposes. A. D. Petford has devised a new photo-electric method of measurement 
(14) which could be used here to great advantage. In fact, the visual and photo- 
electric measurements made on plate E28 of our Table I are in very good agree- 
ment, as his Fig. 3 shows. The photoelectric measurements however show 
considerably less scatter in repeated settings and are indisputably quicker and 
easier to obtain. 

The solar line displacements corrected for instrumental effects, as explained 
above, were then expressed as Doppler shifts, (1 mm displacement = 9:19 km/sec) 
and the results are shown in Fig. 1. Fig. 1 (a) is the mean of the direct and reversed 
readings for the February and April measurements and gives the observed curve 
of the small spot of Plate E28, while Fig. 1 (b) refers to a similar mean of the large 
spot of E22. 

The ordinates in Fig. 1 are arbitrarily set equal to zero for the photosphere 
on the solar centre side of the spot and expressed in km sec~!, while the abscissa 
is given by the spectrum height in mm. The approximate extent of the umbra 
and penumbra is indicated for both spots. 

2.1. Reduction to sunspot velocity fields.—In the reduction of the observed 
velocities, we have followed the scheme proposed by Professor Plaskett (15) and 
followed by T.D. Kinman(10,11). The first step in deducing the motion charac- 
teristic of the spot is to remove the effect of the observer’s motion and to allow for 
the displacements which we know will arise from solar rotation and limb effect. 
The Mercury computor of the Oxford University Computing Laboratory was 
used to calculate the corrections at representative points up the slit height, and 
interpolations were made for the remaining measured points. The corrections for 
the observer’s motion may be calculated immediately from the observed epoch 
of exposure and the polar coordinates of the slit. For solar rotation at latitude B, 
we have used (2) 

R,=2(1—0°229 sin? B)kmsec~. 
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The limb effect corrections (2) were found to be negligible for the disk positions 
considered. ‘These corrections show what velocity shift would occur across the 
slit height in the undisturbed photosphere. The lines in Fig. 1 show the run of 
these velocity corrections across the slit and are positioned for the small and large 
spots to give zero photospheric velocity on either side of the spot. It is the shifts 
above and below these lines which we believe represent the true sight-line velo- 


cities in the spot. 
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To deduce the spatial motion in the spot region from the sight-line velocities, 
we transform the coordinates of the measured points on the projected disk to 
polar coordinates (r, ¢) on the solar surface and origin at the centre of the umbra. 
In this system, r is the horizontal radius from the centre of the umbra and ¢ is the 
position angle measured counter clockwise from the direction of solar rotation. 
If AB and Ad are the differential heliographic coordinates measured from the 
centre of the umbra (heliographic coordinates B, L) then: 


r= R,[(AAcos B)? + (AB)?]*? and tand= AB 


AA cos B 
where Ry = radius of the Sun in kilometres. 

Motion in the sunspot is conveniently considered in terms of the cylindrical 
velocity components u=f?, v=rd, w=, which are the radial, tangential and 
vertical velocity components respectively. If can be shown (11) that the sight- 
line velocity V is given in terms of these components by the linear relation 


V =u(cos cos y, + sing cos y,) + v (cos cos y, — sin P COs y,) + WCOS yz 
where cos y; is a direction cosine which depends upon the position of the points of 
measurement on the disk and is defined (15) as follows: 


cos y, = sin (L — L,) cos By sin (8 + 6,) cosec 8 
COS y, = [sin B cos By cos (L — Ly) — cos Bsin Bo] sin (6 + 6,) cosec 0 
Cos yz = — cos (8+ 6,) 


where 0, ~ Cs S = angular semi-diameter of Sun 
Po 


and sin (0+ 6,)~ 3 


Ideally, we should have liked to make a complete least squares solution for the 
three unknowns u, v, w, but our observational material was insufficient to provide 
a satisfactory set of equations. (It is hoped to attempt such a complete solution 
in future investigations.) We have therefore based our calculations on Kinman’s 
conclusions, which were obtained from 12 least squares solutions for u, v, w, at 
four values of r for each of three plates. These showed that the tangential 
component vw is entirely due to random error and that the vertical component w, 
relative to the photosphere, is zero. 

As a first description of the velocity field, we have then assumed that the motion 
is entirely a radial flow outward from the centre of the umbra (v=o, w=0). 
Fig. 2 (a) shows the results of the calculation for the small spot of E28, and Fig. 2 (b) 
is a similar plot for the large spot of E22. Positive velocities are now shown for 
both sides of the spot, since both refer to outflow velocities. For the small spot, 
where we have four sections of the same spot, we may go further and attempt a 
solution for another component of motion. Assuming that the spot has radial 
symmetry about the centre of the umbra, the two sides of the spot on four plates 
lead to eight observational equations for the state of motion at any one r value. 
Still assuming v = 0 (10, 17), we have used these equations to obtain a least squares 
solution for finding components u and w. These solutions were carried out for 
eight r values and lead to the results shown in Fig. 3. 

We find as Kinman did, that w is zero within the limits of error. We see 
also that the values of the radial velocity thus obtained for the four plates are 
quite similar to those shown in Fig. 3 (a) for plate E28 only. 
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2.2. Errors.—The error in setting on the solar lines after taking direct and dif 
reversed readings and the mean of two sets of readings is approximately 10 microns, sp 
which is equivalent to a velocity error of 0-1 km/sec. The comparison of photo- it 





electric and visual measurements on one plate indicated a velocity difference of 
about 0-15 km/sec in the penumbra (the most difficult place to make a setting) and 
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the discrepancy was considerably less in the region away from the penumbra. 
7 We can thus consider a line setting error 0-15 km/sec throughout as a generous 
estimate. No correction has been applied to the observed sight-line velocities 
for scattered light, but Kinman (10) gives an estimate for a reduction in the 
observed penumbral velocities of about g per cent due to thiseffect. However, as 
Plate 6 (a) shows, the definition of our spectra is good and observations have shown 
that the instrumental scatter for the 35 m telescope is very similar to that for the 
19 m telescope used by Kinman. 
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Discussion.—The general form of motion in the small spot is very like that 
obtained by Kinman, and within the limits of our observational material, it 
again appears that the flow is wholly radial. The maximum radial velocity 
(mean uncorrected value of 1-1 kmsec™) is achieved towards the edge of the 
penumbra, and we see that the flow extends well beyond the penumbral region, 
falling to zero at approximately 30 00o0km from the spot centre. Considering 
Kinman’s linear relation (11) between maximum radial velocity and spot size, 
we see that there is fair agreement for the small spot, but our maximum radial ‘ 
velocity for the large spot (mean uncorrected value approx. o-6kmsec~') falls 
well below the value (approx. 3-5 kmsec~!) which may be estimated from the 
linear relation. In addition, the velocity field for the large spot is of quite a 
and different form from that of athe small spot. The large spot pattern is much more 
pas, spread out, falling to zero at approximately 70 ooo kms from the spot centre, and 
oto- it does not show such well defined maxima in the penumbra. 

e of The high dispersion spectra (0:17 A/mm) reveal a phenomenon (discovered at 
and the McMath-Hulbert Observatory (12)) which is unobservable in the low 
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dispersion spectra. This feature, which appears in Plate 6(b) for the large and 
the small spot, will be referred to as “‘line-flare’’. The appearance of the plate 
indeed suggests that the main line is undisturbed and that the measured shift is 
caused by this diffuse ‘‘flare’’. The line-flare now observed extends over some 
1z20mA from the centre of the undisturbed line, corresponding to a sight-line 
velocity up to 6kmsec~. Evershed (9) first observed that spectral lines showed 
diffusive widening in the direction of displacement near the penumbral limits, but 
he remarked that this phenomenon might well be involved with Zeeman widening. 
The line-flare cannot be due to Zeeman splitting here, but is probably a 
Doppler effect due to turbulence, with a mean motion of matter away from the 
umbra (as shown by the flare appearing on opposite sides of the line on either side 
of the umbra). The same feature has been recently noted in Bumba’s investiga- 
tion (3) of spot velocity distributions. A complete account of his work (4) has 
come to our notice during the preparation of this paper. He also suggests that 
the line-flares (‘‘ flags ’’ in his notation) are a result of Doppler shifts, and that the 
line as a whole is not shifted. In low dispersion work, such as Kinman’s, the 
line-flare and the undisturbed line would be blended and consequently the 
amalgamated flare would contribute in a different way to the measured Doppler 
shift. This phenomenon could well account for the differences observed in 
velocity displacements in high and low dispersion investigations. In measuring 
the high dispersion plates, it is difficult to say exactly to what part of the line the 
micrometer setting refers, but presumably (Plate 6(5)) it is closer to the 
undisturbed line than would be possible with lower dispersion. 

The plates used here for radial velocities are not calibrated for photometry, 
but high dispersion calibrated plates have been obtained in the early summer of this 
year for a full-scale study of the intensity and extent of line-flare, and to investigate 
more fully the part it plays in our measurement of sunspot velocity fields. Further 
possible considerations in such investigations are the effects of spot age and disk 
position (13, 17) on the sunspot velocity components. 
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Summary 


Relative oscillator strengths of 200 lines of Ti1t and Tir have been 
measured from spectra of a vortex-stabilized arc running in titanium tetra- 
chloride at 10080°K. Comparison is made between the present measure- 
ments and those of other workers; these data are summarized on a scale 
which is approximately absolute. It is shown in the case of titanium that the 
excitation potential effect discovered by Allen and Asaad is merely a selection 
effect. 





1. Introduction.—The measurement of oscillator strengths in the laboratory 
continues to be an essential prerequisite to the quantitative interpretation of 
astronomical spectra. The intensity of a spectrum line emitted or absorbed by 
an element in given physical conditions depends not only on N, the number of 
atoms present, but on gfN, where g is the statistical weight of the initial level, and 
f is the oscillator strength of the transition. In general, g is known accurately 
from the level designation, but not f. In complex spectra detailed values of gf 
may be obtained only by measurement. The present paper describes some gf 
measurements in titanium, which is of great astrophysical importance. 

2. Equipment.—The main piece of apparatus used was a modification of the 
“‘water-stabilized arc’’, originally developed by Maecker (13) and Jiirgens (8). 
The basic unit, which we shall call the ‘‘ vortex’’, is drawn diagrammatically in 
Fig. 1. It is constructed from a short solid brass cylinder down the axis of which 
is a hole of diameter 8 mm, widening into a central chamber of diameter 16 mm. 
Liquid is injected tangentially into the chamber through two orifices; there it 
forms a little whirlpool and then flows out through the ends of the axial hole, where 
it spreads, without splashing, over the ends of the vortex, which are shaped like 
telephone earpieces. An 80amp. D.C. arc, powered by a heavy generator, was 
struck down the axis of the vortex between two hollow carbon electrodes, and the 
plasma was observed end-on through one of the electrodes. 

The advantages claimed for this type of arc are three-fold: 

(a) Stability. Although the stability of the vortex arc has been questioned 
(12) there can be little doubt that when it is running correctly it is far more stable 
than an ordinary arc in air. 

(5) Composition. If the arc is run at a known temperature and pressure in 
a liquid of known chemical composition, the composition of the plasma may be 
computed. 

(c) Temperature. Since the plasma is constricted by the vortex, the tem- 
perature is high. Further, since the outside of the plasma is cooled by the liquid, 
thus reducing the electrical and thermal conductivities, most of the current is 
carried in the middle of the plasma and the heat generated there does not escape. 
The middle becomes extremely hot, thus permitting the excitation of high-level 
lines. ‘The high radial temperature gradient is not necessarily an advantage, 
however. 
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Fic. 1.—The vortex-stabilized arc. 


It might be thought that the excitation of a metallic spectrum could be achieved 
simply by running the arc in a solution of salt. Experiment shows that in these 
circumstances the salt does not penetrate suitably to the plasma (12, 24). How- 
ever, one may obtain a metallic spectrum by using a liquid which contains the 
metal in molecules held together by covalent or coordinate bonds. Examples of 
such compounds are the organometallic compounds such as Pb(C,H;), lead 
tetraethyl, the carbonyls such as Fe(CO), iron pentacarbonyl, and some chlorides 
such as TiCl, titanium tetrachloride. Such compounds are usually rare and 
dangerous. Exhaustive search in the chemical literature revealed that titanium 
tetrachloride was almost the only suitable liquid, and it is a happy coincidence 
that titanium is of considerable astrophysical interest. ‘The nearest rival to 
titanium tetrachloride appears to be SnCl, stannic chloride, which is more 
expensive but less objectionable. The tetrachlorides of vanadium and germanium 
might also be suitable, but on the whole it seems that the future of vortex arc lies 
with the non-metals, for which there are many suitable liquid compounds. 

The spectrum was recorded photographically with a Littrow-type quartz- 
prism spectrograph having a dispersion of 16 A/mm at Hy and 6 A/mm at 3200A. 
The plates used were 1oin.x 4in. Ilford Long-Range Spectrographic. Line 
intensities were compared with the spectrum of the standard radiator described 
by Euler (6). This is a carbon arc whose positive pole is at 3995°K and whose 
emissivity has been given by Euler. Measurements were made with a Hilger 
H670 recording microphotometer. In order to calibrate the plates, a rhodium- 
on-quartz step filter, having densities of o to 1-2 dex (i.e. in 10-based logarithms) 
in steps of 0-2 dex, was placed in front of the spectrograph slit and the spectrum 
of the carbon arc was photographed. The lines of the titanium spectrum have an 
intensity range of more than 1-2 dex and, to calibrate the plates over the whole 
intensity range, three exposures of 1, 8 and 32sec were made. The resulting 
partial characteristic curves were combined into one complete curve. 
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3. Experiment.—Preliminary experience with the vortex arc was obtained by 
running it in water and in carbon tetrachloride. Experiments in which salts 
were dissolved in water or in carbon tetrachloride produced only a few feeble 
metallic lines. 

In 1959 June the arc was run for the first time in pure titanium tetrachloride. 
Dense white fumes were produced as soon as the liquid came in contact with the 
air. ‘The vortex arc seemed to run fairly well except that there was quite a lot of 
light outside the vortex. This suggested that some arcing might have taken 
place between an electrode and the vortex instead of from electrode to electrode 
through the vortex. When the apparatus was examined afterwards a quantity of 
a friable yellow solid, probably a mixture of titanium dioxide and titanium oxy- 
chloride, was found to be clogging up the vortex and the space between vortex 
andelectrodes. It was thought that most of this was not present while the spectro- 
graphic exposures were being made, but was formed shortly afterwards when the 
water vapour in the air had had sufficient time to react with the titanium tetra- 
chloride. The spectrum showed many hundreds of neutral and ionized titanium 
lines but no chlorine lines. The titanium lines were rather diffuse and in many 
cases strongly self-absorbed. 

During the course of the next few months many more spectra were taken of 
the vortex running in pure TiCl, On some occasions it was known that 
undesirable arcing between electrodes and vortex had occurred, whilst on other 


ved occasions the arc was thought to have run correctly. The spectra usually had 
ese approximately the appearance described above, but sometimes when the arc 
w= was running incorrectly the spectrum showed a strong continuum crossed by 
the absorption lines. Because of the uncertainty in knowing how the arc would run 
s of with this liquid which hydrolyses so readily when exposed to the air, it was 
ead decided to use a solution of titanium tetrachloride in carbon tetrachloride. ‘The 
des solution which was found to be most suitable was one part of titanium tetra- 
and chloride to four parts of carbon tetrachloride. The tendency for the above- 
um mentioned yellow solid to form and short the arc was considerably reduced. 

nce When the arc was running in such a solution it was possible to approach fairly 
| to closely and examine the plasma. On nearly every occasion the arc was running 
ore satisfactorily. ‘The spectrum lines were much sharper than in the case of pure 
jum TiCl, and were rarely self-reversed. The stronger lines were, however, self- 
lies absorbed, and all lines suffered from the general absorption of the fumes. Lines 


of neutral chlorine were present and they were noticeably broader than the 
rtz- titanium lines. 





0A. 4. Composition of the plasma.—The composition of a plasma resulting from 
sine the evaporation of a mixture of composition aCCl, + TiCl, is given by the solution 
bed of the following seven equations : 
108e , 
Iger m+ Na+ Nz + Mgt M+ Mtn a. (1) 
um- 1 2 3 4 5 6 e kT 
ms) 4(a+1)n,+4(a+1)n, —- M— Me =o (2) 
rum 4(a+1)ng+4(a+1)my— an,— an, =o (3) 
e an No + nN, + n.—n, =o (4) 
hole Sy — NN =o (5) 
Sons — nM, (6) 
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where 7,, Mg, Mg, 14, M5, Mg, are the numbers per unit volume of the atoms Ti, Tit, 
C, C+, Cl, Cl* and n, is the number per unit volume of electrons. P is the total 
pressure (equal to atmospheric pressure in the present work), T the absolute 
temperature and k Boltzmann’s constant. Sx denotes the Saha function for 
the element X. 

Equation (1) is a statement of the perfect gas law. Equations (2) and (3) 
follow from the molecular composition of TiCl, and CCl, Equation (4) 
expresses the fact that the plasma is electrically neutral. Equations (5), (6), (7) 
are Saha’s equations for Ti, C, Cl. 














Fic. 2.—Concentration of neutral and ionized titanium resulting from the evaporation of a 
solution of composition TiCl,: CCl,=1: 4 by volume. 


The Saha function is given by 


log S= —(x—Ay)O— Hog 0 +27-24+ log “+175 metres~* (8) 
0 
where @ = 204° x is the ionization potential in volts, and u,(7) and u)(T) are the 


7 ’ 
partition functions of the ion and neutral atom. Partition functions were 
obtained from Claas (5) for the lower temperatures and by calculation for the 
higher temperatures. Unsdld (25) has shown that the ionization potential 
should be lowered by 

Ay =7 x 107°n,!8 volts 
where n, is in electrons per metre*. Ay was at first neglected in the equations 
(5~7) and an approximate value of n, was obtained. The Saha functions were 
then adjusted and the equations were solved again. A further iteration was 
thought to be unnecessary. 

The numbers per unit volume of neutral and ionized titanium atoms and of 
electrons are shown in Fig. 2 for the mixed plasma. It should be noted that 
equations (2, 3) and Fig. 2 are valid only if the TiCl, and CCl, evaporated into 
the plasma at the same rate. 

5. Analysis of the spectra.—Six spectra were chosen for measurement and 
analysis, all of the mixed TiCl,+CCl, plasma. Absolute intensity measure- 
ments were made of as many Ti and Tit lines, free from blending, as was 
possible. Some of these lines had previously been measured by the Kings (9, 
10), others were new. 
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Measurement of the temperature was one of the most important and difficult 


tasks. 


Three methods were used by the writer. Comparison of his intensities 


with the titanium oscillator strengths published by the Kings (9, 10) gave the 


excitation temperature. 


The ionization temperature was measured by com- 


paring the intensities of Tit and Tim lines. A further check on the order of 
magnitude came from an examination of the variation of intensity of 'Ti1 lines 


from the centre of the edge of the plasma. 


The first of these methods must be 


regarded as the most important; the other two provide good confirmatory 


evidence. 


The intensities of emission lines are related to their oscillator strengths by 


3 
log = =const.—OV 
gf 


(9) 


where A is the wavelength of a line, J its intensity, and V the excitation potential 
In Fig. 3 we have converted 


of the upper level in volts. 


As before, © = 5040/T. 


King’s Ti 11 gf-values to relative intensities at four different temperatures, and 
have compared them with our measured intensities from spectrum No. QP 145(1). 
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Fic. 3.—Comparison of measured log intensities with intensities calculated at different tempera- 


tures from King’s f-values. 


sterad~'; the abscissa gives the logarithm of the calculated intensity in arbitrary units. 


Least scatter is expected to occur at the temperature of our experiment. 


The ordinate gives the logarithm of the measured intensity in watt m-* 


At 


lower temperatures, lines of high excitation potential (E.P.), marked by crosses in 
Fig. 3, will be displaced to the left; at higher temperatures, to the right. We see 
that our temperature cannot have been far from 10080°K. ‘The Tit lines give 


3 
the same temperature. A direct plot of id : V, which, for optically thin lines, 


should be a straight line of slope —@, confirms this result. It should be noted 
that any error in King’s determination of his furnace temperature will be per- 
petuated by a similar error in our temperature. 
that the experiments of the Kings are still the best ever done on titanium oscillator 
strengths ; this does not obscure the fact that independent temperature measure- 


ment is most desirable. 


The writer is of the opinion 


It may be noted that R. B. King (11) still believes 


“that no very large systematic error was made in the furnace temperature measure- 
ment’’, and further evidence for this is presented in the section of this paper 
dealing with the excitation potential effect. 

In order to calculate the ionization temperature it is necessary to put King’s 


gf-values for Tit and Ti tt on a single scale. 


it was possible to put King’s gf-values on the absolute scale. 
cussed in the section on absolute values. 
given by the solution of: 


(07 


of 


Using the formulae of Allen (3), 
This will be dis- 


The ionization temperature is then 


3 
- {log (=) +0V,} =logn,—2724+ 15 logO+OV 
4 


(10) 
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where the suffixes o and + refer to the neutral atom and the ion, V, and V,, are 
the excitation potentials of the upper levels of an atomic and ionic line, V is the 
ionization potential, and m, is the number of electrons per cubic metre. Mean 
values of the braces on the left-hand-side were calculated for all measured atomic 
and ionic lines at various temperatures. The required temperature is obtained 
when the difference between the mean values of the braces is equal to the right- 
hand-side. Ionization temperatures obtained in this way were in satisfactory 
agreement with excitation temperatures for the three spectra for which such an 
analysis was possible. The degree of agreement is indicated by the following 
table: 


Spectrum No. Excitation temperature Ionization temperature 
QP 145(1) 10080 °K g500 °K 

QP 145(2) 9500 9500 

QP 146B 9500 gg0o 


The principle of the third method of temperature measurement was to measure 
the change of intensity of neutral titanium lines from the cool perimeter of the 
plasma to the hot interior. With increasing temperature, the intensities of the 
neutral lines at first increase because of increasing excitation. The intensities 
reach a maximum and then decrease because of increasing ionization. From a 
measurement of the ratio of the maximum intensity to the intensity in the middle 
of the plasma, it is possible to deduce the central temperature. Certain difficulties 
in this measurement made it possible to obtain only a lower limit to the central 
temperature; thus the central temperature of QP 145(1) was found to be at 
least g500 °K. 

6. Results and comparison with previous observations.—Relative log gf-values 
were calculated from measured intensities by equation (9), assuming tempera- 
tures from g500°K to 10080°K for the different plates. The scales on which 
the gf-values were calculated were those of King and King (9) for Tir and ef 
King (10) for Ti 11, but they have been reduced to absolute values as described 
in the next section. Comparison of our values with previously published values 
has been made in detail, and this has been satisfactory. In particular the standard 
deviation of our log gf-values from the Kings’ values is about + 0-2 dex in Til 
and +o-1 dexin Till. 200 lines were measured in all, of which 113 are not on 
the Kings’ lists. 

Oscillator strengths for Tit have been published previously by King and 
King (9), van Stekelenburg (22), Ostrovskii (18), Hefferlin et al, (7), Mitrofanova 
(15), Allen (1) and Boyarchuk and Boyarchuk (4); and for Ti 11 by King (10), 
van Stekelenburg (22), Mitrofanova (15), Rountree (20), Allen (1) and Boyarchuk 
and Boyarchuk (4). The published values are of variable quality and nearly 
all on different scales. In view of this and of the 1955 I.A.U. Report on line 
intensities (14), which stresses the need to ‘‘ make available all theoretical and 
experimental data from various sources critically balanced and combined ’’, the 
writer has compiled a table (23) which lists all published titanium oscillator 
strengths known to him. This table gives 884 measurements of 420 lines in 
137 Til multiplets and 369 measurements of 202 lines in 67 Ti multiplets, 
all reduced to a scale which is approximately absolute, as described in section 7. 
The table gives also the weighted mean of all the measurements. 

The present paper gives only those lines which have been measured by the 
writer. Tables I and II give the results for Ti1 and Tit respectively. ‘The 
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are successive columns are the multiplet number and wavelength written in the same 
the order as in the Revised Multiplet Table (16), our measurement of log gf, and the 

ean weighted mean of all published values, omitting the astrophysical measurements 
mic of Allen (1) and Boyarchuk (4) and the laboratory values of Mitrofanova (15). 
ned 
ght TaBLe I 
tory 
1 an Absol . : 
: solute log gf-values in Tit 
ving 
Mult. A Tatum Mean Mult. A Tatum Mean 
ture 6 4681°9 — 1°32 — 1°09 115 3789°3 —0°65 —0°65 
4667-6 — 1°42 —1'24 3798°3 — 1°02 — 1°02 
9 4112°7 — 1°64 — 1°63 126 4820°4 — 0°47 —0°47 
12 3998°6 —o'12 —0'06 129 4186°1 —0°22 — 0°22 
3989°8 —o'18 —0'20 130 3919'8 — 1°22 — 1°22 
sure 3981°8 — 0°27 — 0°33 131 3724°6 —0°85 —0-08 
the 3964°3 — 1°08 — 114 133 3547°0 —o'80 —o'80 
the 3962-9 —1'18 —I°ls 145 4617°3 +0°25 +0°29 
ities 13 39582 —O'14 —o'18 4623°1 —o-'06 —o'ol 
3956°3 —0°45 — 0°45 4629°3 — 0°44 —0°44 
m a 3948°7 — 0°46 —0'48 4650°0 — 0°94 —0°76 
ddle 3929°9 — 1°00 — 1°07 46452 —0'92 —o'76 
Ities 14 3947°8 —o°81 — 0°93 146 4481°3 — 0°35 +004 
tral 3921°4 — 5°52 —1°§2 4496°1 — 0°23 —0°23 
a 17 3752°9 —0°23 —o'lo 148 4276°4 —0°46 —0°46 
3729°8 — 0°36 —0°39 157 4885°1 + 0°03 +0°22 
37717 —0'74 —0°99 4899°9 — 0°07 +o°10 
lues 18 3689°9 — 1°03 — 1°24 161 4404°9 +0°08 +0°08 
yera- 3669-0 — 1°34 — 1°37 162 42631 +o'18 +018 
hich 19 3653°5 —0'13 +o'11 4282°7 — 0°09 — 0°09 
4 of 3642°7 —o'lg +o°o1 4265°7 —o'80 —o'80 
jis 3635°5 —O1r , —0°O4 185 4030°5 — 0°03 — 0°03 
ibed 3671°7 —o'9I — 1°07 4021°8 —0'50 —o'50 
ilues 24 3354°6 —0°55 —0'13 189 3822°0 —0°25 —0'25 
dard 3341°9 —0'99 — 0°39 206 4151°0 —0'02 — 0°02 
Tit 27. = 3199°9 +0°05 +0°12 4159°6 — 0°43 — 0°43 
wpe 38 4981°7 —oO'l4 +o°61 207 4099'2 — 0°42 — 0°42 
5007°2 —o'41 +0°64 220 4203°5 —0°45 —0°45 
42 = 4555°5 — 0°59 —0'52 4200°8 — 0°65 — 0°65 
and 4544°7 — 0-62 —0'55 4188-7 —0°49 — 0°49 
nova 4512°7 — 0°52 —0°45 221 4154°9 — 102 — 102 
(10), 4518-0 — 0°39 — 0°33 231 4856°0 +0°37 +0°37 
4522°8 —0°46 — 0°37 233 4742'8 +0°16 +o0°12 
chuk 
44  4305°9 +0°35 +0°41 246 3938-0 —0°33 — 0°33 
early 53 4840°9 —0°75 —o°62 252 4256°0 —0'56 —0'56 
line 56 3904°8 +0°13 +0°14 253 4137°3 —0'26 —o'26 
and 57 3786°0 +0°17 +002 4143°1 — 0°42 — 0°42 , 
, the 59 3598°7 — 0°92 — 0°92 4131°2 — 0°69 — 0°69 
lator 62 3292°1 —o'731 —0°23 254 4071°2 — 0°36 — 0°36 
3 75 4698°8 — 1°07 — 1°04 260 4805°4 +0°73 +0°73 
es in 4722°6 —1'07 — 129 4792°5 +0°40 +040 
lets, 80 4078°5 —o'1o —o'16 284 4237°9 0°00 0°00 
on 7. 4082°5 —o'81 —0°74 291 4278°2 +or'rr +Oor'll 
83 3710°0 —o'31 —0°39 292 3926°3 +026 +0°33 
113 4457°4 +0°15 +0°15 296 4127°5 +0°26 +0°48 
y the pais bas ine 
4455°3 0°22 0°22 301 4224°8 —0°28 0°28 


The 






Mult. 


10 


II 
12 


13 
15 


16 
19 


20 


26 


31 


A 
3361-2 
3372'8 
3383°8 
3380°3 
3387°8 
3394°6 
3407°2 
3409°8 
3234°5 
3236°6 
3242°0 
3254°3 
3252°9 
3251°9 
3217°1 
3222°8 
3197°5 
3157°4 
3130°8 
3444°3 
3461°5 
3477°2 
3491! 
3322°9 
3329°5 
3335°2 
3340°3 
3343°8 
3346°7 
3308°8 
3318-0 
3326°8 
3168-5 
3162°6 
3161°8 
3161°2 
4012°4 
3813°4 
38146 
3721°6 
3587°1 
3561°6 
3573°7 
3341°9 
4395°0 
4443°8 
4450°5 
4294'1 
4337°9 
4344°3 
3190°1 
3202°5 
4468°5 


4501°3 


Tatum 


—ovol 
—0°'05 
— 0°26 
— 0°69 
— 0°65 
— 0°74 
— 1°89 
— 1°84 
+0°04 
—o'll 
—0°40 
—0°65 
—0'56 
—o-72 
—0'°52 
—o-50 
— 1°52 
—I°41 
—048 
— 0°93 
— 0°97 
— 1°04 
—1'17 
— 0°39 
—0°48 
—o'58 
—0'72 
—1°'2I1 
—1'14 
—1'21 
—I1'I4 
—1'I2 
— 0°36 
— 0°49 
—o088 
—1'07 
=s39 
Felt 
=a oe 
—o-98 
— 1°43 
—1°81 
— 1°36 
— 0°04 
—0'30 
— 0°64 
— 1°46 
—0'92 
—o'89 
— 1°94 
+oOr°1! 
+o°1o 
—0'52 
—o°79 


F. 


B. Tatum 


TaBLe II 


Absolute log gf-values in Tit 


Mean 
+0°05 
—0°'05 
—o'1g 
—o-71 
— 0°66 
—0°76 
— 1°88 
—1°89 
+0°07 
—0'07 
=—O35 
—0°70 
—o'60 
—o'72 
—0'58 
—0°74 
— 1°52 
—I°41 
—0'48 
—o'9gI 
— 0°99 
— 1°06 
—1'21 
=—O°33 
— 0°48 
— 0-60 
— 0°74 
— 1°24 
— 1°16 
— 1°23 
— 1°16 
— 1°20 
—0°36 
—o'50 
—0'76 
—0'92 
= 5°99 
=F 35 
=i93 
— 1°04 
ao 
—1°81 
— 1°36 
— 0°04 
—0-30 
— 0°64 
— 1°46 
—0'92 
—o'89 
— 1°94 
+o°15 
+o°10 
—0°52 
—0o°79 


Mult. 


32 
34 


36 
40 
41 


45 
50 


51 
52 


72 


73 
75 


82 
84 
87 
88 
89 
94 
98 


99 


105 


107 
114 
115 


125 





A 


4341°4 
3900°5 
3913°5 
393270 
3232°3 
4417°7 
4312°9 
4307°9 
4315°0 
3276°8 
4563°8 
4590°0 
4399°8 
3641°3 
3624°8 
3402°4 
3352°1 
3337°9 
3332°1 
3321°7 
3315°3 
3248°6 
3261°6 
3271°7 
3278°3 
3282°3 
3741°6 
3776°1 
3706°2 
3659°8 
3662-2 
3679°7 
4549°6 
4572°0 
3224°2 
32183 
4028°3 
4053°8 
3504°9 
3510°8 
3287°7 
4350°8 
3535°4 
3520°3 
34564 
3452°5 
3465°6 
4163°6 
4171°9 
3748-0 
4911'2 
4488-3 
44I1I°I 
3480°8 


Tatum 
— 2°22 
—o0'08 
—0'20 
—5 so 
—0'06 
—I1'ls 
—I'I4 
— 1°00 
—0'96 
— 0°68 
—0°84 
—1°81 
— 1°20 
—Or'rs 
—o'16 
— 0-98 
— 1°02 
—1°I2 
— 0°04 
—0°30 
— 0°64 
+0°24 
+0°32 
+0o°oI 
—oO'17 
—0'23 
+0°27 
— 1°04 
—0'38 
—0'26 
—o0'28 
— 1°25 
—O'l4 
— 0°24 
+0°16 
+0°13 
—o-'81 
— 0-98 
+0°59 
+026 
+0°38 
— 1°69 
+0°24 
— 0°02 
— 0°04 
— 0°53 
—0'84 
+0o°oI 
—O'ls 
—0'24 
— 0°42 
—0'46 
— 060 
+ 0°32 





Mean 


— 2°22 
—0'08 
—0'20 
baa 
—o'12 
—I'ls 
—I'l4 
— 1°00 
—0'96 
—0°68 
—0°84 
—1°81 
— 1°20 
—O'15 
—o'16 
—0'98 
— 1°02 
—I1'I2 
— 0°06 
—o-30 
— 0°64 
+0°24 
+0°24 
+o°oI 
—o'18 
—0°25 
+0°27 
— 1°04 
—0-38 
—0'26 
— 0°28 
— 1°25 
—o'14 
—0'24 
+0o'16 
+0°13 
—o'81 
—0'98 
+0°59 
+0°26 
+0°38 
— 1°69 
+0°24 
—0'02 
— 0°04 
— 0°53 
— 0°84 
+0o°oI 
—O'1s 
—0'24 
— 0°42 
—0°46 
— 060 
+0°32 
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7. Absolute gf-values.—It has not been possible to obtain absolute gf-values 
from the present experiment. Line intensities were measured in absolute units 
since the absolute intensity of the carbon arc standard radiator was known. It 
was possible to calculate the numbers per unit volume of atoms and ions of each 
species in the plasma, and since it is believed that most of the lines were optically 
thin it should therefore have been possible to calculate absolute gf-values. When 
this was done, however, values were obtained which were of the order of about 
10~* of what one might reasonably expect. There are two possible reasons for 
this discrepancy. ‘The rate of evaporation of the titanium tetrachloride into the 
plasma may have been much lower than that of the carbon tetrachloride, so that 
the concentration of titanium atoms was much less than calculated. On the 
other hand the very dense fumes which arose from the liquid during the experi- 
ments may have caused serious loss of light. The writer believes that fumes 
accounted for by far the greater part of the light loss, and that differential evapor- 
ation was negligible; but the matter is not yet settled. Comparison of the 
present log gf-values with those of the Kings shows that any wavelength depen- 
dence of the extinction due to the fumes was certainly less than 0-15 dex in the 
range 3000-5000 A, and therefore this caused no serious systematic wavelength 
dependence of the present log gf-values. 

According to Allen (3), the conversion factor from the scale of King and King 
for Tir to absolute values is given by: 


log gf (abs.) = log gf (rel.) — 3-22. 
Allen does not give a conversion factor for Ti 11 but it is possible to calculate one 
from his empirical formulae. The most reliable of these formulae is the one for 


d"s—d"p multiplets, viz: 
Laf=0'2338ut 


where ut is the statistical weight of the upper term. The multiplets 1, 2, 13, 14, 
16, 24, 26, 65, 88, 98, 99 are all d*s—d*p multiplets; none of them is intercom- 
bination or interparental, and the writer has been able to estimate total multiplet 
strengths on King’s relative scale for them, either from his own or from published 
data. The mean difference between the calculated absolute log Xgf and the 
measured relative log Xgf is 2-37+0°15. The small standard deviation gives an 
indication of the reliability of Allen’s formula in this case. 'The conversion from 
King’s relative scale for Ti 11 (10) to the absolute scale is therefore: 


log gf (abs.) =log gf (rel.) — 2°37. 

Absolute gf-values calculated on these bases were used to compute an ioni- 
zation temperature for the vortex (Section 5), and the agreement between this 
and the excitation temperature seems to confirm the correctness of both the 
excitation temperature and Allen’s formula. 

8. Oscillator strengths and excitation potential.—It was pointed out by Allen 
and Asaad (2) that there is a tendency for experimental oscillator strengths of 
iron-group atoms to increase with excitation potential. I have examined the 
effect in titanium. 

There are several things which could cause the excitation potential effect and 
we may group them under the following headings: 

A. The effect may be a systematic experimental error, such as: 

(i) Wrong temperature estimation. 
(ii) Failure to allow for self-absorption. 
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B. It may be a selection effect. This has been suggested by Traving (26). 
C. It may be real in either of the following senses: 

(i) There happens to be a strong array with high E.P. 

(ii) Oscillator strength is a function of E.P. 


Before we consider the effects to which these might give rise, we notice the 
following. Underestimation of temperature would produce an apparent increase 
of gf with E.P.; overestimation would produce a decrease. If A(i) is the explan- 
ation then we must believe that no experimenters have overestimated but all have 
underestimated their temperatures, and indeed very grossly underestimated 
them. A(ii) could give rise to the E.P. effect, if self-absorption, acting only on 
low-lying lines, were neglected. ‘This would affect only emission experiments, 
whereas the E.P. effect is known to occur also conspicuously in absorption experi- 
ments. A selection effect, in which lines of high E.P. and low gf are too faint to 
measure, must operate; the only question is how important selection is. King 
(x1) has pointed out that the array d"s—d"p, which might be expected to have 
strong lines, has many high-lying multiplets, and this undoubtedly contributes 
to the E.P. effect. A similar suggestion has been made by Allen (26). C(i) and 
B are not entirely separable. As for C(ii) there is no theoretical reason for 
believing gf to be a function of E.P. and we must consider this only after we have 
eliminated all other possibilities. ‘Thus we are led at the outset strongly to 
suspect that the E.P. effect is a selection effect. 

We can predict qualitatively and quantitatively what type of E.P. effect A 
and B would give rise to, on the assumption that gf is not a function of E.P. Let 
us suppose that, for all E.P., all gf-values lie roughly between the limits « and f. 
Fig. 4 shows how we might expect measured values of loggf apparently to vary 
with E.P. In Fig. 4 (4) is shown the expected increase of log gf with E.P. if the 
temperature were underestimated. For emission experiments the abscissa is 
the upper E.P., for absorption experiments it is the lower E.P. Fig 4 (a) shows 
how one might underestimate log gf for lines ending at the ground term if one 
did not allow for self-absorption in these lines. This would occur only in an 
emission experiment; the abscissa in Fig. 4 (a) is the lower E.P. In Fig. 4(c), 
in which the abscissa is as in Fig. 4 (b), we see the effect of selection. ‘The upper 
limit of measured oscillator strengths remains the same, but for high E.P.’s we 
cannot measure the lines of small oscillator strength. We can measure lines 
right up to the point of the triangular area in Fig. 4(c), but not beyond. We 
may write equation (9) as: 


log gf = log A°J — const. + OV 


and we see that lines of slope © are lines of constant A°J_ or roughly of constant 
intensity. ‘Thus we see that the hypotenuse of Fig. 4(c) represents the lower 


~ . Fora number of 





limit of measurable intensity and should have a slope of 


reasons the lower limit should not be expected to be as clear-cut as indicated, but 
nevertheless selection should give rise to an E.P. effect roughly as shown. 

Fig. 5 shows the E.P. effect actually obtained by various workers. A line of 
slope © is drawn in each diagram for comparison with the observed distribution 
of log gf-values. These graphs speak for themselves and indicate that selection 
is the main cause of the E.P. effect in titanium. 
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The solar and stellar gf-values of Allen (1) and Boyarchuk (4) require additional 
comment, however. The ratio of the solar and stellar values to laboratory values 
shows a slight residual dependence on E.P. which cannot be explained by 
selection. Thus log gf (Allen) —log gf (laboratory) decreases slightly with E.P. 
Allen used © =org for the reciprocal excitation temperature in the Sun. If we 
adjust his gf-values to @ = 1:0 the residual E.P. dependence disappears. Log gf 
(Boyarchuk) —log gf (laboratory) increases slightly with E.P. Boyarchuk used 
@=1'1 forthe Sun. If we adjust his gf-values to © = 1-0, again the residual E.P. 
dependence disappears. 


























log gf log of log of 
Pr KK 3 wee B | 
ae “ _-_-— “ 
lower E.P. EP. ‘EP. 
(a) (4) (c) 


Fic. 4.—Apparent increase of log gf with E.P., due to: 


(a) Non-allowance for self-absorption in the low-level lines 
(b) Underestimation of temperature 
(c) Observational selection. 
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Fic. 5.—Observed increase of log gf with E.P. in volts. The observers are respectively Tatum 
(this work), King and King (9), and Allen (x). The straight lines are drawn corresponding to the 
temperature quoted by each author. 


The corrections are small, however, and there is no doubt that the major 
reason for the apparent increase of oscillator strength with excitation potential is 
observational selection. This result is important in that it shows that published 
f-values in titanium are not marred by serious systematic error (reference (15) is 
an exception) and are therefore trustworthy. It is expected that the same 
would be found for other atoms in the iron group. 

g. Comparison with theory.—There are two classifications of the spectrum of 
Tit extant; that of Russell (in 17) and that of Rohrlich (19). The main 
differences in these classifications are in the configuration assignments of the 





322 j. B. Tatum Vol. 122 


singlet and triplet terms of odd parity. Rohrlich computed multiplet strengths 
based on his assignments, and it is accordingly of much interest to compare his 
multiplet strengths with the observed values. He computed strengths not only 
for multiplets allowed by the LS-coupling rules, but also for some multiplets 
which are, on his assignment, interparental; he did not include configuration 
interaction and did not compute strengths of intercombination multiplets. Before 
we compare laboratory values with Rohrlich’s values, we should note that the 
purpose for which Rohrlich computed multiplet strengths was primarily to 
guide him in making his assignments and not necessarily to produce accurate 
multiplet strengths for astrophysical purposes ; in his own words, referring more 
particularly to interparental multiplets, ‘‘ these values are only very rough approxi- 
mations ’’. 
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Fic. 6.—Comparison of observed multiplet strengths with the calculated strengths of Rohrlich, 
for three different arrays. The abscissa is the observed log Z gf on the absolute scale. The ordinate 
is Rohrlich’s calculated log X gf on a relative scale which is different for each array. The 
left-hand scale refers to d’s—d*p and d*s*—d*sp; the right-hand scale refers to d*s—d*sp. 

The comparison has been made in Table III and Fig. 6. The multiplets are 
grouped into the three arrays of Rohrlich’s classification. The columns of the 
table are: multiplet number from the R.M.T. (16) and the multiplet designation; 
the logarithm of the multiplet oscillator strength, log 2gf, calculated from 
Rohrlich’s multiplet strengths; the logarithm of the observed multiplet oscillator 
strength from (23). It should be noted that the Rohrlich values are relative 
and are on different scales for each array; the observed values are intended to be 
absolute. A colon after the multiplet designation indicates that Russell’s 
assignment (i.e. Russell’s 1947 assignment given by Moore (17), not the 1927 
Russell assignment (21) which is given in Rohrlich’s paper) differs from 
Rohrlich’s. A comma indicates that the multiplet is interparental on Rohrlich’s 
assignment. A semi-colon indicates both. Note also that in multiplet 284, the 
upper term designation is u1D°, not v!D® as given in the R.M.T. (16). The 
reason has been stated by Moore (17). In Fig. 6, singlets, triplets and quintets 
are indicated by circles, triangles and squares. Underlined multiplets are those 
for which Russell’s assignment differs from Rohrlich’s. Open symbols are inter- 
parental on Rohrlich’s assignment, full symbols are L.S-permitted. 

The trends which appear from Fig. 6 seem to indicate that for L.S-permitted 
multiplets, the agreement between Rohrlich’s values and experimental values 
is tolerably good, including those multiplets for which Rohrlich’s and Russell’s 
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Taste III 
Total multiplet strengths 
Multiplet Rohrlich Laboratory Multiplet Rohrlich Laboratory 
3d°4s—3d°4p 115 ; +r +180 — 0°29 
238 «6s G-z'H®: +109 §=©—0°36 © ore » +134 +026 
ae onan 156 wv F°?, +1°84 +0°26 
240 ~y : 9 3 5F_-5? 42:38 — 1°30 
284 b'D-u1D° +0-75 0°00 37 ary"; 2°3 13 
286 aH-y'G®; —o'52 +0°34 
288 —z1]° +1°40 +0'50 3d74s°-3d74s4p 
291 —x1H° +1°24 +or1r ID_-1N0 a mee 
157 aG-y*H® = +1°74—S ss +-:0°66 . - papi ae Ss aan 
160 —v®G° +1°73 +0°32 53 -y1D°; éine ail 
162 -t°F° +1°61 +0°59 6 1F0 : ; : 
Pie ae 5 ~y +154 +0°14 
220 BP-®P®; +042 +0°13 57 aa ot owen 
231 a®H-2*1° +1:90 +0°77 58 ~yip?: +18 —om 
233 -@H® +1°76 +0°81 59 et 5 42 —o-92 
260 c®P-s*D® +127. +1°04 126 aG-y'!F° eee ~063 
8 OES +1 +8 131 -IG°; +168 +031 
42 ou. +2 +20 4 @F-8F*, +155 —0-44 
44 -D® +175 +0°94 5 a, +0 on 
145 aP-w*D® +171 +0°77 6 3G + 67-07% 
146 7". +t. 6 +eR 12 ~F°; +072 +0°36 
3T)0. . . 
34°4s—3d*4s4p rl re : ; oe ’ a 
244 b'G-y'H®; +0°96 — 0°25 19 y®G®:; +1°84 +0°53 
287 aH-x1G®; +2°04 +0°30 23 —w*D®: +162 +o°14 
289 -—~y'H®; +1°90 +0°52 24 —-®G®; | +1°31 +041 
102 BF-x°F®; +2°03 —0'57 26 —~v'D®; +061 — 1°60 
104 -y®G®; +0°90 +0°13 69 a®P-2z°S° +o°71 —0'97 
106 —uw®D®; +1°86 +018 72 -~+'°D°®; +048 —0°46 
108 -x8G®; +0°86 — 086 75 —w®D®; +0°47 —046 
109 vw'D®; +1°95 —o'lo 83 -y°S°, +0°89 +0°06 


assignments differ. An exception is multiplet 59, which has, however, low 
experimental weight. The scatter of points is worse for the interparental 
multiplets, although there is no great systematic difference from interparental 
to permitted. Many of the interparental multiplets in Fig. 6 have different 
assignments on Rohrlich’s and Russell’s classifications. ‘This makes it difficult 
to decide whether the large scatter is due to the ‘‘ very rough approximations ”’ 
of Rohrlich’s interparental multiplet strengths or to inappropriate assignments. 
The most serious discrepancies are for multiplets 37 and 57, both of which have 
high experimental weights. ‘The reason for the discrepancy in multiplet 37 is 
not clear. Multiplet 57 is a1D—z1P®. Rohrlich assigns z 1P° to 3d74s(b?P)4p 
and therefore the multiplet is interparental and should have zero intensity. 
Russell assigns z 1P° to 3d24s(a?D)4p; the multiplet would therefore be permitted 
and strong, as is observed. It would seem that the latter assignment is more 
appropriate. However, we should then have to follow Russell and assign 
y 1P° to 3d4s(b?P)4p; multiplet 58 would then be interparental, whereas it is 
observed to be strong. This situation is typical of any attempt to improve the 
observed-calculated agreement by changing assignments. It would probably be 
fair to conclude that there is not necessarily a unique ‘‘correct’’ assignment for 
some multiplets. 
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Summary 


Spectral types and luminosities on the revised Yerkes (MK) system are 
given for 6 stars of type O, 201 of type B, and 13 of type A, in the region of the 
Scorpio—Centaurus Association. Microphotometer tracings were used to 
estimate the line ratios. A mean modulus is determined for the certain 
members of the Association, and variations of the mean modulus with type 
and luminosity are discussed briefly. 





As a part of an extensive study of the Scorpio—Centaurus Association, the 
spectral classification of early-type stars in the region of the Association was 
undertaken. The distance modulus obtained from the spectral type and lumino- 
sity is a sensitive criterion for membership of the Association. 

Spectrograms of the brighter stars were obtained with a three-prism spectro- 
graph at the Cassegrain focus of the 30-inch reflector (dispersion 36 A/mm at Hy), 
while for the fainter stars a two-prism spectrograph was used at the Newtonian 
focus of the 74-inch reflector (dispersion go A/mm at Hy). These will be referred 
to as ‘Cassegrain’ and ‘Newtonian’ spectra respectively. The Cassegrain 
spectra were obtained primarily for radial velocity measurements and were too 
narrow for reliable classification by visual inspection, whereas the Newtonian 
spectra were sufficiently wide and well resolved to be classified visually. However, 
in order that the classification system should be as homogeneous as possible, all 
spectra were treated in the same way. ‘They were recorded with a Hilger micro- 
photometer at enlargements of x 7 for Cassegrain and x 15 for Newtonian spectra, 
and relative line intensities were estimated from the tracings. In addition, the 
wings of the hydrogen lines and any peculiar features were examined with a bino- 
cular microscope. ‘To retain homogeneity the classification of Cassegrain spectro- 
grams has been suspended since the recent installation of a new set of camera and 
collimator optics, and will not be resumed until a new library of standard spectra 
has been built up. 

The standard stars for classification were selected from the work of H. L. 
Johnson and W. W. Morgan (1), M. L. Woods (2), and A. de Vaucouleurs (3). 
A comprehensive list of standard stars suitable for use in the southern hemisphere 
has been compiled by W. Buscombe (4). The distribution of the standards 
observed is given in Table I, which shows that the main gaps occur in luminosity 
class II. 
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TaBLe I[ 


Classification standards recorded 


Cassegrain Newtonian 
L I II III IV V I II III IV V 
i 
O7 I 
Oo I I I I 
O9°5 2 2 I 
Bo 2 I I 
Bors I I I 2 I 
Bi 3 2 2 I 
B2 I I 2 I 2 
B25 I I 
B3 I I I I I I I 
Bs 3 I I I I I I 
B6 I I I I 
B7 2 I I I 
B8 I I 2 3 I 2 
Bo I I I 4 I I 2 
Bo's I I 
Ao 2 I 2 2 2 
Al I 2 
A2 I 2 2 I I 
A2'5 I 
A3 I I 
A4 I I I 
AS I I I I I I I I 
AZ I I 
Ag I 
Totals: O 5 Oo 5 
B 45 B 28 
A 22 A 14 
72 47 


The 220 southern stars which have been classified on the MK system (1) 
are listed in Tables II and III. Table II gives those stars for which Harvard 
photovisual magnitudes (5) are available, while the magnitudes of stars in Table 
III have been derived from various other sources. The columns show: 

(1) the star name, or Gould number (6), or, if these are lacking, the 
Durchmusterung number; 

(2) the Henry Draper Catalogue number (7); 

(3) and (4) the galactic longitude and latitude to the nearest hundredth 
of a degree, on the revised system (8) ; 

(5) the spectral type and luminosity class ; 

(6) the Harvard photovisual magnitude, or, in Table III, magnitude derived 
from other sources as shown by reference numbers in the final column; an 
asterisk indicates that the magnitude has been corrected for duplicity by + 0-3, and 
a dagger that it has been corrected by +0°6; 
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(7) the colour excess derived from various sources, all converted to the B-V 
system; the values of the intrinsic colours for MK types were taken from (9), 
(10) and (11): 
(8) the modulus m)—-M=m-—3E,—M, where M is the visual absolute 
magnitude corresponding to the spectral type and luminosity as given by H. L. 
V Johnson and B. Iriarte (12); 
(9) the spectrograph used; C indicates Cassegrain and N Newtonian ; 
I (10) remarks, and references to notes at the end of the Tables. 
I 
I 
Tas.e II 
Star oo + ole 6S m(HPy) E, m—MSp. Notes 
I Lum. ” oa0m - 
72G Col 41534 23843 —23°37 B2IV 5°62 +o:02 886 C 
83G Col 43071 243°55 —23°08 B3V 6°83 “Sn 853 N 
. 9oG Col. 44506 241°63 —20°78 B1Vn 5°86* +0:03 9:37. C 
42G CMa 46189 23617. —16°47 B4IV 5°83 000 68803 CC 
50G CMa 46036 241:'29 —17'54 B&V 5'91* +005 566 C 
: 70G CMa 49028 240:12 —14°67 B8IV 6°43 ooo 7:13 N 
m 30G Pup 49336 247715 —17:19 B3V 6-06 +o05 761 C 
8G CMa 50012 237°51 —12°36 B3IV 6-90 —o'0oz2 940 N 
P 82G CMa 50093 23611 —11°62 B3V 6°24 —-oor 7:94 C 
21° 1655 51200 23315 — 899 B3I1V 7-r1* oi 961 N 
I 
1044G CMa 51925 23812 —1067 B3V 6°19 —o02 789 N 
21° 1695 52437 233°73 — 804 B4Vne 6°33 —o:04 7°73 N Em. HB 
I 27° 3540 52812 23850 —10:°02 B3V 6:66 +003 827 N 
I 26° 3880 54224 23851 — 868 BrV 6°38 +0706 980 C 
124G CMa 54669 23633 — 715 B3V 6°85* -—oor 855 N 
136G CMa 55857 23981 — 7°65 B3IV 5°86 si 8-36 C 
79G Pup 56733 250°09 —11'°92 Bs IV 5°71 +o02 745 C 
80G Pup 56779 24852 —11'12 B2V 5°32* +005 767 C 
1499G CMa 56876 239°73 — 658 Bs5 Vn 634 —oor 7:34 C 
157G CMa 57573 230955 — 412 B3V 6°45 +o:04 8-01 C 
96G Pup 59550 245°05 — 6°52 Ba-s5I1V 6-10* o'00 «6gtcoo =—C 
n (1) y Pup 59635 25164 — 991 B3IV 5°70* +0°04 808 C 
evard 5G Pup 60606 249°85 — 7:96 B3:Vnek 613f +008 7:59 N (a) 
Tabl 18G Pup 60753 262°80 —1442 B6IV 6°70 Lo 8:20 N 
apie 125G Pup 61068 235°53 + 060 Bail 5°66 000 10°56 C 
» the F 26° 4723 61687 24213 — 237 B7IV 6-74 —o0z 7:94 N 
1446G Pup 61925 251°53 — 7°46 B3IV 6:26* +014 834 N 
163G Pup 62753 25433 — 813 B2Vne 6:54 +012 868 C (8) 
redth 55G Car 62758 270°92 —16°71 BsIV 6°63* +005 828 N 
168G Pup 63028 240°55 + 0721 B3IV 6°55 od gos C 
. Cc 
rived N 
1; an Cc 
3, and N 
Cc 


rN) 
aN 






79G 
254G 
258G 

10G 

36° 


° 


35 
87G 


36G 
95G 


50G 
55G 
60G 
61G 
62G 


68G 
| 
82G 
92G 
103G 


124G 
125G 


157G 
166G 


172G 
167G 


186G 
206G 
205G 
238G 
240G 
261G 
264G 
36G 
42G 
54G 
70G 
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Star 


Car 
Pup 
Pup 
Vel 
4359 


4413 
Car 
5245 
Vel 
Car 


Vel 
Vel 
Vel 
Vel 
Vel 


Vel 
Vel 
Car 
Cha 
Vel 


Vel 
Vel 
Car 
Vel 
Cha 


Vel 
Vel 
Vel 
Vel 
Car 


Vel 

Car 
Car 
Cha 
Car 


Car 
Car 
Car 
Mus 
Cen 


Cen 
Cen 
Cen 
Cen 
Cen 


67536 
67698 
67888 
68324 
69106 


69620 
70839 
7IOI5 
72350 
73390 


74067 
74196 
74455 
74535 
74560 


74753 
75149 
75311 
75416 
76161 


77320 
79186 
80094 
80781 
83979 


84816 
85980 
86440 
86466 
86659 


88206 
91533 
91619 
93237 
93607 


93737 
96919 
97583 
99264 
100198 


100673 
100841 
101189 
102776 
102878 


ju 


276711 
242°66 
254°52 
263°33 
254°50 


254°15 
272°87 
252°32 
262°69 
274°15 


260°18 
270°31 
266°59 
279°57 
270°59 


268-10 
265°32 
273°90 
292°40 
267°88 


264°81 
267°36 
277°65 
275°75 
295°58 


271°85 
273°25 
279°34 
278°19 
288°39 


279°13 
285°72 
285°55 
297°19 
289°96 


288-04 
291°14 
292°43 
296°31 
293°45 


291°78 
294°47 
294°38 
296°17 
295°98 
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TaBLe II1—(continued) 


pl 


— 16°14 
+ 4°04 
— 2°64 
= Pee 
= Sa 


— 053 
—I1I'gI 
+ 2°62 
— 3°21 
—10°51 


+ og! 
— 6°85 
— 3°61 
— 6°67 
— 6°67 


— 4°48 
— 1°69 
— 8-42 
— 21°65 
— 2°44 


+ 1°95 
+ 2°24 
— 6°56 
= oa 
— 21°05 


+ 6°66 
+ Pag 
+ o'10 
+ 1°65 
— 11°30 


+ 3°39 
~ 0-50 
— 0°05 
— 18°39 


= 320 
— 10°50 
+ or1O 


+ 6°93 
= 2a 
— o16 
— 1°72 
— o60 


Sp. & 
Lum. 
B4 Vn 
Bs Ve 
Bs III 
B3 V 
Bors III 


B6 V 

B3 III 
B3 Ill 
Bs IV 
B3 Vn 


Bg V 

B8 IV 
B3 Vn 
Bg III 
B4 IV 


Bo Vn 
B2 II 

B2 Vn 
Bo IV 
B6: Vn 


Bs Ia 
B3 IV 


Ao Ia 
Bo Ia 


B3 Ill 
Ao Ia 


Bo V 
Bog II 
Bo IV 
B3 Vne 
Az Ia 


m(HPy) 


6°30 
6°64 
6°30 
5°58* 
7°22* 


6°92 
6°25* 
7°93 
6°37 
5°58* 


5°12 
5°91* 
5°67* 
5°59 
4°92 


ss7" 


5°97 


6°24* 
6°32 
5°42* 


5°92* 
5°69 
3°58 
6°03 
6:22 


5°28* 
6°13 
6°14 
6°36* 
4°98 


6:06 
5°34 
5°40 
5°54 
6°30 


5*00* 
3°24 
5°24 
4°70* 
5°59 


E 


v 


+0°02 


+0°08 
— 0°03 
+0°17 


+0°13 
+0:08 
+o°10 


+0°04 


+0°13 
— 0°03 
0°00 


+0:'02 


+0°06 
+0°45 
+0°'09 
+0°04 
— 0°02 


+0°03 
+0°33 
+0°03 
+ 0°03 


+0°03 
+0°06 
+0°03 
+0°06 
+0°06 


+008 
+026 
+0°52 
+0°13 
+0°04 


+o'18 
+0°24 
+0:06 
+0°21 


0°00 
+o°o1 
+0°09 
+0°04 
+0°23 


m—-MS 


7°64 
7°64 
8-26 
7:28 


11°61 


7°23 
8-91 
9°63 
8-18 
7°16 


4°13 
6°61 
7°37 
5°99 
7°06 


9°59 
gol 
7°04 
5°91 
6-99 


7°98 
9°81 
Fae 
7°43 
7°42 


7°53 
6°91 
7°49 
7°55 
8-24 


7°54 
12°35 
11°58 
6°97 
7°36 


12°52 
11°42 
4°62 
7°81 
13°30 


4°40 
6-21 
5°17 
6:28 
12°10 
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Em. Hf 





118 


132 
39 
41 
42 
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TaBLe II—(continued) 


Sp. & 


I 
HD b Lum. 


mHPy) E, m—MSp. Notes 
104174 —15°62 BoVn 5°25" +011 4°32 
104841 -— 076 B3IV 5°09* +010 7°29 
104878 — 5°83 AoV 5:04 +o1r 4°61 
105382 +11°63 B6III 4°68 —0°03 6°58 
105416 +1357 ArV 5°80* +0°03 4:21 


105509 +1789 AzIIl 5°88 +0°35 3°63 
106490 + 3:78 Baill 2°93 000) «653 
107696 + 4:99 B8Vp 5°78* —o-or 5°68 
108114 +27°38 BglII 5°75 +0°03 6:06 
109668 — 629 B2IV 3°09* +0°03 6°30 


a 
y 
V 


110073 +22°84 B8IV 5°02* +001 5°69 
110335 + 317 B7IV +o'10 5°82 
110461 + 6:90 BoV ° +oO'II 5°50 
110506 + 667 BoVn , +0°04 5°42 
110879 - B2'5 V +0°05 5°37 


111774 ‘ B7 V ? +0°07 6:26 
111904 , Bo Ia +0°36 11°46 
112409 : B8 V . +012 5°02 
113120 Br Vne +0°22 9°05 
113314 Ao V +014 3°82 


-_~ 
® 
~ 


113703 Bs IV +oor 6°88 
114529 B8 V +oor 4°82 
114911 B8 V ° 0°00 65°04 
115823 Bs III ° +0°03 8:26 
116072 B4 Vn : +0°04 7°94 


116226 B7 IV ° +0°02 7°47 
118978 ; Bo IV +0°04 5°81 
118991 B8 V +012 5°14 
IIQI59 B2IV +0°08 9°21 
120307 B2IV +0O°OI 7°03 


120642 Bo Vn +0°07 4°83 
121190 B8 V +0°04 5°81 
122879 Bo II , +0°30 

124197 B6 V an 7°58 
125238 B3 IV +o-or 6°49 


125721 B3 V +o-or 7°78 
125823 Bs III — 0°04 
126341 B3 III : +0°02 
126769 B8 IV 5° +0°02 5°88 
126981 B6IV+Ar: +o1r 7:18 





126983 AoV+B con 5°06 
129116 B3 IV +0:02 6:80 
129954 B2V +012 
131058 , B4V +0°05 
131120 B6 V 000 = 580 
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Tas_e [I—(continued) 


Star HD ne pit a m(HPy) E, 
36G Lup 131562 32127. + 5°56 A2V 5°80* 
37G Lup 131657 323°62 + 9:91 BoV 5°98* +009 
40G Cir 135160 31967 — 2°83 BrV 5°80 +0°09 
8 Cir 135240 31969 — 2:90 O8n 5°39* +020 
45G Cir 135591 32012 — 2:64 Oogll 5°66* +0°13 
7G TrA 135737 316°50 — 862 B3V 6°65* +0°08 
e Lup 136504 329°22 +1032 B3IV 4°25T 0*00 
x! Aps 137387 313°85 —13°99 B3IV 5°82* +0°03 
«x? Aps 138800 314°33 —14°42 B8IV 5°65 +0°05 
2G Nor 139129 326°98 + 2°52 BoV 5°40 +018 
94G Lib 139160 343°25 +2331 BolIV 6°30* +0°05 
12G Nor 141168 327°33 + 080 B8&V 5°86 +008 
19G TrA 142139 323°84 — 5°33 ArV 6:21* +0°04 
6G Sco 142250 345°55 +20°01 B7V 601 +008 
7G Sco 142301 347712 +21°51 B8IV 5°85 0°00 
116G Lib 142378 351°62 +25°66 Bs V 6:20* +0°14 
10G Sco 142883 350°90 +24:09 B3V 5°88 +0°20 
11G Sco 142990 34812 +21:21 B8&V 5°65* —0-04 
7 Lup 143118 338°76 +11°03 B2V 3°86* +001 
8 Sco 143275 350°09 +22°50 BoV Z311f +0°14 
23G TrA 143448 324°54 — 5°95 B3IV 6°84 +0°04 
151G Lup 143699 33912 +1043 BsIV 5°26* +0-o1 
B? Sco 144218 353°18 +23°61 B3V 4°98 +0°14 
6 Lup 1442904 340°82 +11°33 B3IVn 4°27 +o°o1 
19G Sco 144334 350°35 +2085 B8&V 5°87 —o'ol 
w Nor 144480 326°79 — 4°53 BoV 5°71 +0°12 
24G Sco 144661 350°00 +19°98 B7V 6°31 +003 
v? Sco 145501 354°61 +22°70 AolV 6°12 +0'04 
v' Sco 145502 354°61 +22°70 B3IV 4°50* +0°22 
41G Nor 145782 327746 — 511 AolV 6"10* +0°09 
41G Sco 145792 350°99 +19°04 B6V 6°31 +o'18 
46G Sco 146001 g50°38 +1813 B8IV 6:09 +016 
52G Nor 147152 334°01 + 0113 B6IV 5°36 xe 
a Sco 147165 351°31 +17°01 Brill 3°63* +0°33 
66G Sco 147628 34287 + 827 B8IV 5°69* +0°04 
p? Oph 147933 353°68 +1769 B31V 5°42* +0°35 
p? Oph 147934 353°67 +1768 B3V 6°12* ae 
55G Nor 147977 32810 — 677 BglII 5°98* +013 
59G Nor 148379 337'23 + 1°57 Bala 5°61* +0°84 
p Nor 149038 339°38 + 2°51 Bors Iab 5°38* +026 
1G Ara 149485 32723 — 9°35 B8&V 6°11 +0°04 
7 TrA 149671 321°70 —14:28 B7IV 591 +0°07 
62G Nor 149711 340°38 + 2°36 B3IV 5°99 +o'18 
97G Sco 151804 343°59 + 1°95 O8fk 5°59* +0°26 


41° 1021 152235 343°31 + 110 Borsla 6°03 +0°78 





mo>—-M Sp. Notes 


3°90 
5S‘11 
9°13 
9°99 
11°37 


8-11 
6°75 
8-23 
6:20 
4°26 


6°35 
5°52 
4°59 
6°17 
6°55 


6°78 
6°98 
5'S5 
6°33 
7°09 


9°22 
7°93 
6:26 
6°74 
5°77 


4°75 
6°62 
5°70 
6°34 
5°53 


6°47 
6°31 
6°86 
7°04 
6°27 


6°87 
7°82 
5°99 
10°19 
10°90 


5°89 
6-90 
7°95 
Io°’o!I 
10°49 
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47G 
53G 


179G 
191G 
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Star 


47G 
53G 


179G 
191G 


Sco 
Ara 
Sco 
14091 
Sco 


Sco 
Oph 
13072 
Sgr 
4478 


4484 
Sgr 
13804 
Ser 
13816 


4533 
Ara 
Ser 
Sgr 
Ser 


4736 
Ser 
4896 
14462 
5310 


12565 
Sgr 
5005 
CrA 
Tel 


Ser 
Sgr 


Star 


58° 
18G 


969 
Vel 
4212 
1513 
4267 


3017 
3038 
3193 
4080 
1777 


HD 


152408 
153716 
156325 
158186 
158320 


158408 
159975 
160124 
161756 
163892 


164002 
164.402 
164492 
164637 
164816 


164833 
165024 
165402 
165516 


165793 


166286 
167771 
168021 
170978 
173375 


173994 
175191 
176630 
178628 
180183 


182180 
183133 


HD 


62805 
69168 
71216 
71634 
72555 


96248 
96446 
97557 
105580 
111558 


ju 


343°55 
331°88 
353°75 
355°90 
35418 


351°27 
17°00 
356°43 
1°99 
7°15 


712 
7°14 
7°00 
7°34 
6°05 


7°34 
343°32 
20°29 
8-92 
355°74 


13°42 
12°69 
12°70 

9°37 
16°59 


348°45 
9°56 
28°54 
358°54 
340°88 


10°68 
23°35 


jit 


270°96 
262°54 
25881 
273°32 
264°83 


289°92 
290°08 
290°76 
297°66 
30291 


TaBLe II—(continued) 


pli 


1°06 
9°86 
2°93 
1°61 
0°25 


t+++1+ 


+ 12°35 
— 0°67 
+ 048 


+ 0°47 
— oo! 
— 0°23 
— 022 


— 046 
— 13°80 
+ 6°14 
— 0°43 


~ 522g 
— 3°12 
— 1°44 


; —19°61 


— 12°42 
— oe 
— 20°23 
— 25°93 


— 19°03 
— 14°90 


Sp. & 
Lum. 
O7f 

Bs V 

B6 IV 
B3 IV 
B4 IV 


Bo Ia 

Bors II 
B8 IV 
Bir Ibk 
Br Ibk 


Bo II 
O8k 
Bo Ibk 
B3 IV 
B6 V 


B8 Vk 
B4IV 
B4 IV 
B8 III 
B3 IV 


Bs IV 
Bs V 


TaB_e III 


Sp. & 
Lum. 
Bog V 
B3 IV 
B8 V 
B7 IV 
B4V 


Br II 
B2 IIIp 
B3 IV 
B6 V 
B7 Ia 
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m(HP,) 


5°88 
5°74 
6°34 
6°81 
6°67 


3°04* 
4°53 
7°03 
6°46* 
7°44" 


7'45* 
5°72 
6°91 
6:89* 
7°34 


7°16* 
3°76 
6:06* 
6:26 
6°50 


7°58 
6-67* 
6°68* 
6°75 
7°06 


737 
2°07 
7°69 
6-48* 
6°83 


6-20* 
6-99* 


Ey 


+0°38 
+0°08 
ost 
+ol3 


—0°02 
+0°09 
+0°04 
+ 0°32 
+0°40 


+0°23 
+012 
+0°40 
+o°12 
+0°22 


+0°14 
+0°14 


+o°21 
+0°05 


+052 
+o°9r 
+0°32 
+018 
+0°26 


—0°02 
+0°25 
+o°12 


+0°04 


+034 
+0°09 
+o:16 
+0°14 
+or1l 


mo—-M Sp. Notes 


10°24 
6°50 
6-91 
8-92 
8-87 


8-94 
4°66 
9°41 
8-00 
11°24 


11°36 
11°56 
I1‘21 
11°83 
11°08 


12°34 
8-84 
6°76 
11°73 
12°45 


11°72 
10°94 
II‘92 
8-71 
6-98 


7°07 
4°27 
9°14 
7°12 
9°33 


7°88 
7°99 


10°8 
10°0 
9°2 
73 
13°6 
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TaBLe [1 ]—(continued) 


Sp. & 

Star HD [fl = FS my E, mo—M Sp. Notes 
66° 2171 115846 305°80 — 4:83 B4V 71 ete 85 N (7) 
73° «1184 119109 306°70 —11'15 B7V 70 wae 7°4 N (7) 
65° 2652 124182 311'21 — 462 B6IV 7:0 eee 85 N (7) 
16G Cir 128293 312°72 — 7°42 Bs5:V:e 670 +016 7:2 1C(xs)(f) 
55° 6509 136003 322°68 + o92 Balil 6°72 +031 94 N (x3) (14) 
41° 10925 I51SIS 342°81 + 1°73 B3IlI 74 ane 11°8 N (7) 
22° 4500 164359 770 «6+. 036 ~=6=BollIk 8-2 +0°30 12°7 N (7) 
18° 4871 167411 12°70 — 068 Bollk 8-6 +047 129 N (9) 

Notes 


(a) HD 60606. Emission in Hf, Hy and A 4267. 

(6) HD 62753. Emission in Hf, Hy, HS and He. Very broad, badly-defined lines. 

(c) HD 67698. Emission in Hf and Hy. 

(d) HD 107696. A 4072 is extremely strong for the type, and appears the same on four 
plates of this star. 

(e) HD 113120. Emission in Hf and Hy. 

(f) HD 128293. Emission in Hf, Hy, HS and He. Both class and luminosity are 


uncertain. ‘There may also be emission at A 3933, but this region is not well-resolved on 
the spectra available. 


The majority of the colours were found in the work of Schilt and Jackson (16) 
and the Cape Mimeograms (17). Some were taken from the work of Oosterhoff 
(13), Hoffleit (14), Stoy (15), Hogg (18), Hiltner (19), Stebbins et al. (20), and 
Cousins (21). All conversions to the B—V system were linear, with the exception 
of the Schilt and Jackson values, which were converted by comparison with colours 
derived from other sources. A recent paper by T. Walraven and J. H. Walraven 
(22) and a private communication from R. H. Hardie and D. L. Crawford arrived 
too late to be conveniently incorporated. However, it was noted in both cases 
that the colours given were in good agreement with those from other sources. 

The criteria used for type and luminosity classification were the same as those 
quoted by Mrs de Vaucouleurs (3). For B stars the most prominent criterion of 
spectral type is the rapid decrease in strength, from Bz2 to Bg, of the Het line at 
4026. An estimate of the strength of this line compared with that of Hé quickly 
gives a good first approximation for the type. 

Some stars from Mrs de Vaucouleurs’s list have been re-classified as spectra 
of improved quality were obtained. In spite of the poor quality of some spectra 
used by Mrs de Vaucouleurs, there is little change in the assigned type and lumino- 
sity. A new mean modulus for the stars with normal spectra which were con- 
sidered to be certain members of the Scorpio—Centaurus Association by Blaauw (23) 
has been computed using some revised classifications and the improved absolute 
magnitude values (12). The moduli for the individual stars are given in Table IV. 
Three stars, HD 105435, 120324, and 151890, which have peculiar spectra, have 
been omitted, but none of the stars with normal spectra has a distance modulus 
which deviates excessively from the mean on Chauvenet’s criterion. 
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TaBLe IV 


Certain members of the Scorpio—Centaurus Association 


HD Sp. m,—-M HD Sp. myo—-M 
79351 B2 IV 7°18 130807 Bé III 6°33 
85980 B4V 6-91 132200 B2 III 6°79 
103884 B3 V 7°39 132955 B4 IV 7°39 
105382 Bé III 6°58 136298 B3 IV 6°13 
106983 B3 IV 6°55 136664 B3 IV 7°34 
109668 B2 IV 6°30 138764 B6 IV 711 
112078 Bs Vn 6:00 140008 B6 V 6°08 
113703 Bs IV 6°88 142669 B3 IV 6°47 
120307 B2IV 7°03 143118 B2V 6°33 
120709 Bs IV 6°30 147165 Br III 7°04 
120908 Bs V 6°37 148703 B2V 6°93 
120955 Bs IV 6°74 149438 Bo V 7°15 
121743 B2V 6°73 151985 B2 IV 6°88 
122980 B2V 6°88 153716 Bs V 6°50 
125823 Bs Ill 6°71 172910 B3 V 6°46 
129056 Bx III 7°34 175362 B8 IV 6:12 
129116 B3 IV 6°80 


Mean m)—M=6-72 +005 p.e. 


Division into luminosity classes and intervals of spectral type leads to the 
following mean moduli: 


III IV. V Bo-2 B3-4 Bs-8 

m,—M 6-80 6°75 6°64 6°85 6°83 6°50 
p.e., + 0°13 0°07 0:08 0:06 o'10 o'10 
n 6 15 12 12 9 12 


This shows that the use of the new calibration of absolute magnitudes for the MK 
types has eliminated the systematic luminosity effect noted by Mrs de Vaucouleurs. 
The barely significant excess in the mean distance modulus for the B stars of 
earliest sub-classes may indicate that the adopted intrinsic colours are not suffi- 
ciently blue. On the other hand, Blaauw (24) has suggested that the stars 
embedded in the interstellar clouds of Scorpius and Ophiuchus have evolved more 
toward the right from the main sequence in the Hertzsprung—Russell diagram. 
Two of the stars earlier than Bz, HD 129056 and HD 147165, are pulsating stars 
of the BCMatype. ‘There is as yet no secure evidence that the distance from 
the Sun varies significantly from its mean value of 220 parsecs in any section of 
the 60° of galactic longitude subtended by the Association. 

From the distance moduli it will be seen that many of the early-type stars 
examined in this paper lie beyond the Scorpio—Centaurus Association, while some 
of the late B and early A stars are in the foreground. Further studies of the radial 
velocities of possible member stars are in progress. 
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ON TWO METHODS OF BROWN AND SHOOK 
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Summary 


A method given by Brown and Shook for the elimination of short- 
period terms from the disturbing function is essentially an extension of one 
given by von Zeipel. A slight further extension is given here. It is shown 
that the method of Brouwer for the motion of artificial satellites is a case of 
that of Brown and Shook. 

Another method of Brown and Shook uses the orbital true longitude as 
independent variable. ‘The motion of an artificial satellite appeared to 
provide a particularly suitable case for the application of this method, but on 
trial it proved to be extremely laborious. 

Some suggestions are offered with regard to sign conventions in 
Hamiltonian theory. 





1.—Brouwer’s solution (1959) of the problem of the motion of an artificial 
satellite uses Hamilton—Delaunay theory, and eliminates terms of short period 
by a method of von Zeipel (1916). In the application concerned this method is 
equivalent to one given by E. W. Brown and C. A. Shook (1933), which however 
appears to be more general and to need less preliminary explanation. I shall 
restate it here in a slightly more general form than that given by Brown and 


958; 








Shook. 
Let the Hamiltonian be, 
H=H,(p)+m>,H,(p) cos (/,.9,+ nt), (1) 
with 
dp_ _oH dq_ 0H 
a" "3 2° (2) 
and msmall. Neglect of m gives the first approximation 
dq, pe 0H, a dp, _ 
‘dt sag Op, mi a,(p), dt =e (3) 
Write 
oe a,+n,= Vs (4) 


and consider the transforming function 
, , H, / > 
J= Pr Ur — m>. aH?) sin (1-59 + n,t) (5) 
summation being over any set of terms with v,#0, and usually also not with v, 
small. Then 





p, = ad /2g, =p,’ — mE, AP?) 1, cos (i,49, + n.t) (6) 
Vs 


Uy , oH, é be 
ay = I /py! = q — mE FP)“ sin (lg + mt (0) 


r 8 
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and the new Hamiltonian is 


H! = H+ 2 = Hyp) + mSH,(P) 08 (lade + A) 
— m¥/H,(p')=*c0s (,.9, + mt). (8) 


Substitution for p in H, and using (6) gives 


Ho(p) = Ho(p’) — Symay Sy A291, ,c0s (iraq, + nt) + O(m*). (9) 


Then in H’ the terms under the sign >’ of order mare 
m’H,(p') cos (1,4, + nt) — m¥’H,(p) = cos (I,,4, + mt) 
Vs 


= m> i. ay y SAP), ,c0s (1% + n,t). (10) 


On summation with regard to k and reference to (4) we see that all these terms 
cancel except for the difference between p, and p’,. Hence 
H’ = H,(p’) + XomH,(p') cos (/,.9,’ + mt) + O(m*) (11) 

where >, denotes summation over terms in > not included in >’. Thus all 
short-period terms of order m can be removed in one operation and the further 
interest is usually in secular terms and terms of long period. For the latter it is 
well known that a further differentiation with regard to ¢ simplifies the further 
analysis and covers the case of libration. 

If Hy is a function of p, only, so that a,=o for r# 1, and if H does not contain t 
explicitly, so that m,=0, v,is simply /,,.a,, and 


J = py’ q ~ m¥'H,(P') aa (12) 
= Pr’ Gr — = 3’ | H.(0') 08 (L591) 44q- (13) 


In the problem of an artificial satellite, where the Earth’s field is supposed to have 
axial symmetry, q, is the mean longitude and_J is found by direct integration. The 
series can be summed exactly, as is done by Brouwer, who integrates as in (13). 
His transforming function is therefore identical with that found by this method. 
I have verified this directly. 

2.—Brown and Shook give another method, which uses v, the total angular 
distance described, as independent variable. They reduce the problem to a set 
of differential equations, the left sides of which, except for that for t, are all linear 
with constant coefficients in the unknowns, while the right sides all contain 
derivatives of the disturbing function as factors. wu, the reciprocal of the radius 
vector, is treated as one of the unknowns, so that its variation includes those of the 
eccentricity and pericentre. The method has not, so far as I know, been much 
used in planetary theory, apparently on account of the difficulty of expressing the 
coordinates of the disturbing body as harmonic series in v. But it seemed likely 
that the method would be very suitable for the artificial satellite problem, since 
the disturbing function depends only on the position of the satellite itself. The 
method introduced by King-Hele (1958), which gave the short-period terms to the 
first order and the secular motion of the node to the second order, depends on the 
same principles. The method has the further apparent advantage in this problem 
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eliminates H, from the Hamiltonian. 





On two methods of Brown and Shook 337 


that u appears on the right side only in positive powers, except in the equation for t, 
so that the other variables, to any order in J of the theory of the figure of the 
Earth, consist of finite sums of harmonic terms whose coefficients are powers of e. 

When I tried the method I found that it gave the first order terms very easily. 
The second order terms, however, appeared prohibitively difficult. It might be 
expected that, since the expressions on the right sides of the equations all contain 
the factor J, it would be easy to pick out the terms of order J* that are secular or of 
long period. In the equation for u this would imply a search for terms of speed 
near nm to get the motion of the perigee. Actually, the first approximations con- 
tain several terms and the terms on the right are products of three or four such 
sums. ‘To pick out the relevant terms therefore needed all multiplications but the 
last to be carried out infull. The energy equation and an equation given by Brown 
and Shook ((10), p. 174) for the sine of the declination could be used as checks, but 
only at the very end, and would be of little help towards tracing any mistake. 

In comparison with the transformation method, which also involves much 
algebra, the latter has the advantages that at each stage the disturbing function is 
worked out and is used in every equation, and that elimination of short period 
terms simplifies it. In the true longitude method the different partial derivatives 
have to be worked out separately, to the required order, before multiplication by 
the other factors, and become more complicated at each approximation. My 
final conclusicn was that the labour of the method would not be justified by the 
results. 

After this disappointing conclusion I noticed a passage in Brouwer’s obituary 
notice of E. W. Brown (1939), which seems to imply that Brown himself had 
abandoned the method for very similar reasons. As the artificial satellites should 
provide an exceptionally favourable case for the method my attempt may be 
regarded as reinforcing Brown’s conclusion. 

These remarks of course do not apply to the later work of King-Hele and his 
collaborators on atmospheric resistance, where only a first-order theory is involved. 

3.—The standard works on Hamilton—Delaunay theory suffer from a hopeless 
confusion of signs, and it seems worth while to point out that there is a choice 
with, apparently, substantial advantage over the others. In the first place some 
authors write 
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r= Bo (1) 
but almost as many use a negative sign. If the initial coordinates q,. are identified 
with the «, the above gives 8.= —p,9. This is not, however, an argument for 
taking the negative sign. The f, in planetary theory are angles and best regarded 
as analogous to coordinates, and therefore the «, as analogous to momenta. Then 
if a modified Hamiltonian is 


H=H,+mH, (2) 


and Sy, (q, «, t) is the solution of the Hamiltonian equation with H = H), the trans- 
forming function according to the rules 


J=S,0J/da=B, oJ /dq=p, (3) 
— 
H'=H+ > (4) 





So far there is no negative sign. 
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Further, most theoretical writers write further transforming functions as 
J(q, 7, t), leading to the transforming equations 
oJ od 


p _ oq ? p _ aq’ : (5) 
This has an unfortunate consequence. The function J =q,q,’ leads to 
Pr= 4's Pr’ = (6) 


so that coordinate-like and momentum-like variables are interchanged at every 
transformation. To take the transforming function as 


J(q p’, t) or S(q; P, t) (7) 
is however perfectly satisfactory; and 
Jaap, poS, == (8) 
Pp 
leads to gq,’ =4,, p,=p,’, so that the identical transformation does not interchange 
coordinates and momenta. 

One consequence of this form of the theory is that in the equations of motion, 
when the mean longitude / is taken as one coordinate, the principal term in the 
Hamiltonian will be —p?/2Z*, and the disturbing function R, which is an 
addition to the gravitational potential, appears with a negative sign. We could 
compensate this, if we want to do so, by reversing the sign in the definition of H, 
and something could be said for this (it is done in part in Routh’s modified 
Lagrangian function); in works on celestial mechanics, however, it is usually 
achieved by interchanging the p’s and q’s at some stage, not generally with much 
clarity. For myself I doubt whether either method of compensation is worth 
while. 
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THE EFFECT OF TIDAL FRICTION ON ECCENTRICITY AND 
INCLINATION 


Harold Jeffreys 


(Received 1961 January 26) 


Summary 


Darwin’s analysis for the effects of tidal friction on the eccentricity and 
inclination of a satellite’s orbit is adapted to modern laws of imperfect 
elasticity and of friction in ocean currents. It appears that as the conditions 
make the lags in the tides nearly proportional to the speeds his results for 
small viscosity are qualitatively correct, since these have the same property. In 
particular the eccentricity and inclination of the Moon’s orbit should now be 
increasing. 





In a classical paper Darwin (1880, 1908 (esp. p. 342)) studied the variations of 
the orbital elements of the Earth and Moon, and also considered incidentally 
the case of a planet attended by several satellites. His work was quite general, 
but the only case worked out to quantitative results, except for the variation of the 
mean distance and the rotation of the planet, was for a fluid planet of small vis- 
cosity. Since it appears that most of the tidal friction in the Earth-Moon system 
is in turbulent sea currents it has long appeared desirable that this case also should 
be considered ; and since there are now plausible laws of imperfection of elasticity 
in a solid planet they also need examination. 

In the first place we neglect the inclination. The tidal potential due to a body 
at distance R in the plane of the equator is 


U=!OE (f sin® 8.cos2($—0) + § (J — feos" )} (1) 


where v is the true longitude of the body and (r, 0, 6) are polar coordinates. ‘To 
refer to the rotating planet we must replace ¢ by wt +A, where w is the speed of 
rotation and A the longitude of the place. We take 


R=c(1-—acospt), v=nt+fsinft, (2) 
where in elliptic motion, to the first order, 
p=n, «=e, B=2e. (3) 


Then 


2 
U= i [ # sin? {cos (2A + 2wt — 2nt) + ($a + 8) cos (2A+ 2wt — 2nt — pt) 


ce 


+ (3a —B) cos (2A+ 2wt — ant + pt)}+ 4(4 — 3 cos? 0)(1 + 3a cos pt)]. (4) 
The elastic tide raised by a potential 
U=k,r’S, (5) 
in a homogeneous planet of radius a is AbS,, where 


(: $ = 5 pat hy (6) 


19 
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We shall however write it at present by replacing }fma?/Ac* by H, and inserting 
lags 2€, 2€, 2€, 2€, in the periodic terms. If there is imperfect elasticity the 
coefficients in the amplitudes will also need corrections, but these can be adjusted 
later. 

A surface inequality HS, produces an external potential 


=. 
3 {Ma —— (7) 
As the results will be applied to a satellite as ie by the tides raised by another 
I shall not necessarily identify the disturbed satellite with the tide-raising one, and 
consequently state the results in a more general form than is needed for the 
immediate problem. The part due to friction is 


v= — [2 sin? 0 {e) sin (2A+ 2wt — 2nt) + €, (a+ B) sin (2A+ 2wt — ant — pt) 

+ €, (3a — 8) sin (2A+ 2wt — ant + pt)} + $(4 — $ cos? A)e,% sin pt], (8) 
and if accents now refer to the tide-raising body, R and v to the disturbed one, and 
0= 437, 





U,= ae [ey sin (2v — 2n't) + €, (3a’ + PB’) sin (2v—2n't —p’t) 
+ €5( $a’ — B’) sin (2v — 2n’t + pt) +e, a’ sinp’t]. (9) 
For the disturbed body we take 
R=c{1-—ecos(/—w)}, v=/1+2esin (l—w), (10) 


where / is the mean longitude. Rates of change of the elements can be calculated 
as usual. We have 
n*c3 = f{M, (11) 
definingm. Then 
6fMaH 
= 3; 8 
+ €,($a’ + B’ {sin (2/—2n't—p't) + $e sin (3/—m—2n't—p’t) 
— fesin (1+m—2n't—p’'t)} 
+ €, (3«’ — B’ {sin (21—2n't + p’t) + ge sin (3/—mw—2n't + p’t) 
— fesin (l+m—2n't+p’t)} 
+ $e,ex'{sin (I—w + p’t) —sin (l—w—p’t)}]. (12) 
The perturbations are, to the first order, 


ou) (13) 


[e, {sin (2/—2n't) + gesin (3/—mw —2n't) — fesin (l+@— 2n't)} 





i) 


a7 ~ yeoman tat + Fe 
377 I ou, 
dt 4/(fMc) de ~ 


3% 





For a single satellite perturbed by its own tides 
a’=e’, B’=2e’, (14) 
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where we must put e’ =e, c’ =c, 1=nt, p't=1—wm after differentiation. 
Then the secular terms are 
na da 





—< n, Pani 

de 24(f{M)'*aH 

ba = als > (15) 
de___—‘6(fM)'" aH _ 49 

—° ; os e(« Bettatia)| 


The change in dc/dt gives the usual secular variation of mean distance. 
In the case where all speeds are such that the lag is proportional to the speed, 
5 de/dt has the sign of 


~e{(2w — 2n) — ° (2w — 3n) + }(2w—n)— 3n}=e(22w—36n). (16) 


Thus e increases or decreases according as 7 is less or greater than 11w/18. ThiS 
is Darwin’s result. All direct satellites except Phobos should be receding; all 
except Deimos should have increasing eccentricities. 

For the type of linear elastic afterworking that I have considered (Jeffreys 
1957 a, b), the rigidity u is replaced by the operator 
1+1/7'p 
1+1/7p’ (17) 
where p is the Heaviside symbol and 7, 7’ are two constants with the dimensions of 
time, and 7’>7. For a motion with speed o, p must be replaced by ico. In the 
disturbance produced by a tidal potential the elastic properties enter through a 
factor A as in (6), and H contains 





P= Fo 





7’ (tor +1) (18) 
(19449 + 2gpa) T7'to + 19p9 + 2gar7" © 
The lag is proportional to the imaginary part of this. If o7 is large (and therefore 
or’ large) 
19H (7 — 7) 
o(19Hy+ 2gpa)rr" (19) 





The lags are as 1/o. 
For o7 and o7’ small the lag is 
aia 1Qu,o7(1 — 7/7’) (20) 
1QHo7/7’ + 2gpa 
The amplitude is as for rigidity 1) 7/7’ and the lags are proportionaltoc. In both 
cases the lag is small, and the correction to H for imperfect elasticity is O (e*). 

In Darwin’s case the amplitude H is proportional to cos 2¢ and therefore the 
frictional effect to sin 4, and is also proportional to o for small lags. Thus his 
result on the sign of de/dt stands without alteration for the case of o7, a7’ small. 
This however is clearly wrong for the Earth-Moon system since 27/c is of the 
order of 12 or 274 and 7 might be of the order of 204 or 2004. 

In the opposite case where o7 and o7’ are large for all the relevant tides, the 
corresponding factor in de/dt is 


44w* — 63wn + 20n* 3 
(a 3). (21) 
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This has a complicated behaviour on account of the infinities at n= $w, w and 2w. 
Near these the corresponding speeds are small and the approximation fails, but 
it appears that if wr and m7 are both large de/dt must change sign three times as 
n/w varies. But it is worth while to notice that the coefficient is negative for both 
Phobos and Deimos. It is negative whenever n/w <o-2 and therefore for most 
satellites. 

I no longer think that either of the above cases is realized in the solar system, 
but the results may be worth recording. 

For friction in turbulent ocean currents, where one tidal component has greater 
amplitude than the others, it is known (Jeffreys 1959, Section 8.03) that if the 


current is of the form ; 
ven u=asinpt+avsinrpt (v<1) (22) 


the corresponding terms in the frictional resistance following a quadratic law are, 
to order v, 


tka? ( sin pt + 4v sin ret) (23) 


at least as long as rv is small. The relevant terms in the dissipation are in the 
ratio $v*. 

Now it is likely that in the present state of the Earth-Moon system the semi- 
diurnal tides are somewhere near resonance, their amplitudes being largely 
controlled by friction. Then for the ocean €o, €1, €g will be large angles, and since 
the monthly tide produces velocities proportional to its speed it is likely that ¢, will 
contain an extra factor n/w. In that case the tide of speed 2w — 3m will remain 
the dominant one in the effect on the eccentricity, on account of the large 
coefficient of ¢«,, and e should at present be increasing. The complications of 
the problem are great, but the main conclusion seems highly probable. 

For the modified Lomnitz law, in which the creep under constant stress 
increases like ¢*, with « about 0-25 (Jeffreys and Crampin 1961) the lag in a 
periodic disturbance increases by a factor of about 3 over a range of period from 
025% to 100%, ‘This is insufficient to upset the dominance of the term in ¢«. 
Even if 2w—3 is small, this is a case of an exceptionally large lag and the relative 
magnitude of the term is accentuated. 

The general conclusion is that de/dt is in general positive except when n/w 1s 
more than about 2/3. 

This conflicts with the results of G. W. Groves (1960) who finds that e always 
increases. As however he gets de/dc to be independent of w it appears that his 
hypotheses must have oversimplified the problem, since de/dt and dc/dt involve n 
and w in quite different ways. 

The treatment of the inclination is much more complicated, but the essential 
results are in Darwin’s papers. ‘The semidiurnal, diurnal, and fortnightly tides 
all contribute ; the semi-annual tide is negligible. However we see from Darwin’s 
(35), p. 58, that the fortnightly term enters di/dt with a factor sin* i, and what 
matters is the ratio of the lags for the semidiurnal and diurnal tides. For both the 
reasonable laws the ratio of the lags will remain substantially the same as on 
Darwin’s theory. Thus his main result, that the inclination should now be 
increasing, will stand: for small obliquities the obliquity increases if the month is 
more than twice the day, otherwise decreases. 


160 Huntingdon Road, 
Cambridge: 
1961 January 25. 














ut 


th 


ost 


re, 


23) 
the 


ni- 
ely 
ice 
vill 
ain 
rge 


ess 
na 
om 


ive 
» 1S 


ays 
his 


tial 
des 
in’s 
hat 
the 

on 
be 


h is 








No. 4, 1961 Tidal friction on eccentricity and inclination 343 


References 


Darwin, G. H., 1880, Phil. Trans., 171, 713-891. 

—— 1908, Sci. Papers 2, Ch. 6. 

H. Jeffreys, 1957 a, M.N., 117, 506-515. 

1957 b, M.N., 117, 585-589. 

1959, The Earth, C.U.P. 

H. Jeffreys and S. Crampin, 1961, M.N., 121, 571-577. 
G. W. Groves, 1960, M.N., 121, 497-502. 








25 











an 
of 

$0! 
rec 


He 
rai 


th 





DISSIPATIVE INTERACTION BETWEEN SATELLITES 


Harold fFeffreys 


(Received 1961 January 26) 


Summary 


Mutual influence of satellites through tidal friction is considered, and 
with modern evidence on the rate of tidal friction it appears that the effect is 
insufficient to explain approximate commensurabilities of mean motions. 
There is less difficulty in supposing that mutual influence through a resisting 
medium can have been important. 

Kuiper’s theory of the origin of Trojan planets appears to be inconsistent 
with the result that the triangular position is secularly unstable. 





The problem of the effect of the tides raised by one satellite on the motion of 
another was discussed by Darwin (1908, Chapter 8). Since there are several pairs 
of satellites whose mean motions are nearly commensurable, we might expect that 
some dissipative action has produced this state. Any such explanation would 
require the state to be secularly stable. Tidal friction and a resisting medium 
both need to be considered. 

Darwin’s remarks on tidal friction in this connection are rather discouraging. 
He points out that on account of incommensurability of mean motions the tides 
raised by one satellite cannot ordinarily have any secular effect on another; thus 
the orbital elements of all vary independently. ‘The exception would be when the 
mean motions are commensurable, but Darwin says (p. 408) “‘ it is not credible but 
that such relationship should be broken down by tidal friction’’. That is to say, 
he assumes secular instability ; for if the state is secularly stable it will be preserved 
through further changes of the mean distances. 

On both points his arguments need amplification. In the first place, it appears 
impossible to say without detailed examination whether the state of commen- 
surable mean motions is secularly stable or not. In the second, perturbations of 
satellite A by satellite B will affect the tides raised by A, and some of the subsidiary 
tides will have such speeds that they can affect B. It is for this reason that in the 
previous paper (Jeffreys 1961) p was not immediately identified with n. If the 
mean motions of A and B are nm and n’ the mean longitude of A will contain terms 
with p= +2 (n—n’), arising from perturbations, and terms in sin (2/—2n't + pt) in 
(12) of the foregoing paper will not contain n’, and these tide components can 
affect B. What is true, as indeed Darwin points out, is that the amplitude of such 
tides will contain the product of the masses as a factor, and their effect will in 
general be much less than the effect on either satellite of its own tides. Thus the 
orbital elements of the satellites will indeed vary nearly, but not exactly, indepen- 
dently. But this statement need not hold when the mean motions are nearly 
commensurable, since the perturbations of A by B may then be much larger than 
the ratio of the masses of B and the primary. For instance, the mass of Titan is 
only about 1/4000 of that of Saturn, but it forces an eccentricity of about 0-1 in the 


25* 
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orbit of Hyperion. ‘This type of state can occur if the mean motions are such 
that pn — gn’ is small compared with n and n’, where p and g are integers satisfying 
p=q or p=qiti. The former case is that of the Trojan planets. The latter is 
that of the pairs J1, Jir; Ju, Jim; Mimas, Tethys; Enceladus, Dione; Titan, 
Hyperion. 

I have pointed out (1957) that the only satellites whose orbits can have been 
affected to any important extent by tidal friction are the Moon and possibly J1. 
It appears to follow that the commensurabilities, except possibly for the pair 
J1, Ju, cannot have arisen through tidal friction, and a detailed treatment of the 
secular stability seems unnecessary, though it might be relevant to an explanation 
of how the commensurability, once established, could have persisted. 

The effect of a resisting medium depends greatly on the motion of the 
medium. For the solar system as a whole, if there is anything in the nature of a 
pressure, and if the medium is at rest, the density must increase rapidly inwards, 
to such an extent that if the mass within 0-1 astronomical unit of the Sun was of 
the order of the mass of the Sun, that outside the orbit of Mars would be a 
negligible fraction of a gram (Jeffreys 1929, p. 51). I have considered a medium 
effectively in rigid body rotation (1917), but this depends on Jeans’s theorem 
on statistically steady motions, which I maintain to be totally erroneous (Jeffreys 
1937). It would obviously imply that if the rate of rotation was adapted to the 
motion of one planet, the matter at 2° times the distance would have a parabolic 
velocity; and closer to the Sun it meets the same difficulty as a medium at 
rest. Ifa medium is to have an appreciable affect on more than one planet or 
satellite the pressure must be everywhere small, and the steady motion of most of 
the medium must be such that the matter at any point has nearly the velocity of 
a particle in a circular orbit through the point. Near a planet the velocity must 
correspond to that in a variational orbit. 

Since any satellite orbit is essentially an oscillation about a variational orbit it 
appears that secular effects of a resisting medium will be confined to the eccentri- 
cities and inclinations; any effect on the mean distances can arise only through 
the small departures of the medium from the velocity in the variational orbit. 

Interaction between satellites may be due either to viscous or turbulent drag, 
the formation of shock waves, or accretion. ‘The resistance will depend mainly on 
the surface of the bodies, so that it may for small bodies be more important than 
any interaction due to tidal friction. If satellites A, B have different mean motions 
n, n', it appears that,resistance to a perturbation may well give rise to secular terms 
in the mean distances; these must be small, but not so small as for tidal friction, 
for two reasons. First, we have no upper bound to the total amount of the 
effect of a resisting medium, and the fact that most near satellites have small 
eccentricities and inclinations suggests that it has been large. Secondly, the 
fact that surfaces rather than masses will occur in the. action favours it for 
small bodies. 

I considered the effect of a resisting medium on a Trojan planet in three early 
papers (1915a, 1916, 1917). In (1915 a) I took the medium to be at rest. In 
(1917) I took it to be rotating like a rigid body with angular velocity equal to the 
mean motion of Jupiter. In (1916) I took it to have the velocity appropriate to a 
circular orbit. On all hypotheses the amplitude of the faster oscillation about the 
equilateral triangle position would decrease. On those of (1915 a) and (1917) that 
of the slower oscillation would also decrease; on that of (1916) it would increase. 
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But (19154) and (1917) are quite out of the question: (1916) is right, and 
Lagrange’s particles are secularly unstable. 

Now Kuiper (1956) has argued that the primitive planets were much more 
massive than the present ones, and have lost mass; and that Jupiter in consequence 
failed to control some of its outer satellites, which escaped and became Trojans. 
This would require that if a satellite enters the region where Trojan orbits are 
possible, the subsequent oscillations will diminish. But this is contrary to the 
result that the equilateral triangle position is secularly unstable. Even if the 
motion near the third point is largely controlled by Jupiter, there will be some 
influence of the rest of the medium, which would tend to give a compromise 
between the variational motion and that due to the Sun alone. 

Again, if m is the mass of Jupiter, the speed of the slower oscillation is O (m"?). 
Many cases have been worked out where a motion satisfies a differential equation 
of the form 

x+m x(t )x=0 
where x varies slowly (Jeffreys 1915 b; Jeffreys and Jeffreys 1956, Section 17. 13) 
and the general rule is that the speed varies like x? and the amplitude like y~"*. 
That is, if Jupiter has lost mass the amplitude should have increased. 

In fact Kuiper’s hypothesis leads to just the opposite conclusion to what he 
draws. Ifthe existence of the Trojans is to be attributed to change of the mass of 
Jupiter, that change must have been an increase sufficient to counteract the effect 
of a resisting medium. An escaped satellite would become an ordinary asteroid 
and not a Trojan. 
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AN EXAMINATION OF THE STEADY-STATE MODEL IN THE 
LIGHT OF SOME RECENT OBSERVATIONS OF RADIO 
SOURCES 


M. Ryle and R. W. Clarke 
(Received 1961 January 30) 


Summary 


The results of some recent observations of radio sources have been used to 
investigate the “ luminosity function’”’ for the sources. Without making 
any assumptions about their nature, it can be shown from radio observations 
alone that most of the sources lying in any given range of flux density are extra- 
galactic with an emission at 178 Mc/s which exceeds 10*4 watts (c/s)~! ster. 
If it is assumed either that the physical dimensions or that the optical luminosi- 
ties of the sources are comparable with those of the Galaxy, then the emission 
must lie in the range 3 x 107°—10?’ watts (c/s)! ster. 

From these figures it is possible to derive the expected number-flux 
density relationship according to different cosmological models and special 
consideration has been given to the predictions of the steady-state model. 

With the new Cambridge interferometer it has become possible to observe 
sources considerably weaker than those reached in earlier surveys, and hence 
to make a more accurate determination of the actual number-flux density 
relationship; the new data has also allowed more detailed corrections for the 
effects of extended sources and source clustering to be made. 

A comparison of these observational results with those predicted by the 
steady-state model shows a marked discrepancy, the number of sources 
observed with a flux density in the range 0-5 to 2 x 10~** watts (c/s)-' m~* 
being at least 3 +0°5 times that predicted by the model. 

No attempt has been made to select an alternative model to account for the 
observations, but the results appear to provide conclusive evidence against 
the steady-state model. 





1. Introduction.—Previous observations of radio sources situated more than 
20° from the galactic equator have shown there to be an isotropic distribution, 
with some evidence for non-uniformity in depth (Shakeshaft, Ryle, Baldwin, 
Elsmore and Thomson 1954; Mills, Slee and Hill 1958; Edge, Scheuer and 
Shakeshaft 1958; Edge, Shakeshaft, McAdam, Baldwin and Archer 1959). ‘The 
possibility that this latter result is due to effects of the red-shift was first discussed 
by Ryle and Scheuer (1955). 

A new series of observations has now been made at the Mullard Radio 
Astronomy Observatory with an instrument having considerably greater resolving 
power and sensitivity than those used in the previous surveys. 

The results obtained are reliable over a considerably greater range of flux 
density than those of the earlier surveys, and they also allow an examination to be 
made of selection effects which may have been present previously. In particular 
analyses have been made of the distribution of the sources in relation to their 
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surface brightness (Leslie 1961) to find the effect of partial resolution of extended 
sources, and also the effects which might be caused by source clustering (Leslie 
1961 b). 

The combination of the results obtained from ‘‘ whole-sky ’’ surveys made with 
intermediate resolving power (Elsmore, Ryle and Leslie 1959; Leslie 1961) 
with those of a number of limited areas of sky made with the full resolving power of 
the instrument (Scott, Ryle and Hewish 1961) make it possible to obtain a greatly 
improved determination of the number-flux density (N—S) relationship (Scott 
and Ryle 1961). This relationship has been extended to still smaller flux densities 
by applying to the survey of high resolution the statistical method of analysis 
introduced by Scheuer (1957). This technique has been developed by Hewish 
(1961) who has shown that if, as in the present observations, the noise level is small, 
the method can provide information on the N — S relationship even when there is 
about one source per beam area. A number of authors have shown that a study of 
individual sources can only give a reliable estimate of source counts for flux 
densities at which there is not more than one source per 20 beam areas. 

In the present observations it has been possible to extend the analysis to flux 
densities at which there are about 2 x 10* sources per steradian. The results 
confirm the previous conclusion of an apparent non-uniformity in depth. 

In addition to these observations at a frequency of 178 Mc/s, further measure- 
ments of the spectral distribution both of radio sources (Whitfield 1957; Harris 
and Roberts 1960) and of the background emission (Costain 1960) are now avail- 
able which are relevant to the present investigation. 

Before applying the new data to the cosmological problem it is necessary to 
re-examine the question of the radio luminosity function of the sources. Ina 
previous paper (Ryle 1958) it was shown that, even in the absence of optical identi- 
fications, a number of different radio observations could be combined to set 
restrictions on the ranges of luminosity within which the majority of the observed 
sources must lie. Although the number of optical identifications has since been 
increased, it still represents only a small fraction of the sources known; as will be 
shown in Section 2, however, the optical data are valuable in extending the rigorous, 
though more conservative, restrictions on the radio luminosity function set by the 
new radio data. It is now clear that few of the sources are within the Galaxy, and 
that most of those observed are extra-galactic objects having a radio luminosity 
approaching that of the Cygnus-A source. 

The importance of the red-shift in modifying the number-—flux density counts 
is then evident and consideration is given in Section 3 to the predicted curves for 
two particular cosmological models: (a) the steady-state model and (b) the 
Einstein-de Sitter model. 

The results of these computations are compared in Section 4 with the observed 
counts of sources paying particular attention to the observational test of the steady- 
state model. 

2. The derivation of the radio luminosity.—Considerable progress has been 
made in the identification of radio sources with optical objects (Dewhirst 1958, 
1961; Bolton 1960), although the number of identified sources is still only a 
small fraction of the total. It is, however, now clear that most sources are 
extremely faint optical objects (Minkowski 1960; Long and Marks 1961) and 
the observations are consistent with a radio luminosity function in which most of 
the sources observed in a given range of flux density produce a radio emission which 
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is at least 5 per cent that of Cygnus-A. This interpretation is supported by theore- 
tical arguments concerning the mechanism and supply of energy for the emission, 
the gravitational potential of the source, and the minimum optical emission 
expected (Tunmer 1960). 

As will be seen later, the luminosity function derived from optical work is 
entirely consistent with the limits derived from the radio data, but since it must be 
based on a small fraction of the sources the distribution might be seriously in error 
if a class of radio source exists having a low intrinsic optical luminosity. 

The observation of additional sources of small angular size in the area of sky 
covered by the ‘‘spur’’ of galactic emission at /"=30° was first reported by 
Edge et al. (1959) and this result has been confirmed by recent unpublished 
observations at Cambridge. The presence at high galactic latitudes of radio 
sources which are presumably in the Galaxy and which are apparently unrelated 
to optical objects must raise the question of whether there are other optically 
invisible galactic sources. 

A further difficulty is illustrated by the problem of accounting for the source 
Hercules-A, which recent observations have shown to be a double source with an 
angular separation of about 2’ arc between the components (Boishot 1960; 
Williams, Dewhirst and Leslie 1961); the sources may represent a similar 
situation to the Cygnus-A object, although this explanation would require a 
physical separation of the two emitting regions about 200 kpc to either side of the 
optical galaxy. If the intensities of the two components had differed by more than 
3:1 it would have been difficult to recognize the source as double and the apparent 
radio position would have been displaced from the optical galaxy by several times 
the probable error. In such circumstances a very low limit would have been 
wrongly attributed to the optical luminosity of the related object. 

Similar difficulties also apply to the derivation of the distance and hence the 
radio luminosity from measurements of the angular diameter of a source. The 
complex nature and large extent of Hercules-A (and toa lesser degree Cygnus-A) 
suggest that the distances of similar but more remote objects might be seriously 
underestimated. On the other hand a source of small physical size but having a 
large volume emissivity might be interpreted as a more powerful and distant 
source if it were assumed to be of galactic dimensions. 

For these reasons it is important to re-examine the limits which may be set to 
the radio luminosity function from radio observations alone. Having established 
definite but perhaps less stringent limits in this way, the evidence provided by 
optical observations and from theories of the radio emission may allow closer 
restrictions to be set. 


(a) The evidence from radio observations 
In an analysis based on the results of earlier surveys it was shown (Ryle 1958) 
that limits could be set to the radio luminosity (P) from a correlation of 
the following information : 
(i) The number of sources observed in a given range of flux density. 
(ii) The average volume emissivity of different regions of space required to 
account for the observed background emission. 


(iii) The observed isotropic distribution of sources situated at galactic 
latitudes > 20°. 
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Using the earlier observations at a frequency of 81-5 Mc/s, it was shown that (c 
most of the sources must belong either to: 
(i) a class of galactic object having P~ 10" watt (c/s)~? ster! (at 81-5 Mc/s) ar 
and an average space density throughout the galactic halo of p= 10~’ pe~*; or ist 
(ii) a class of powerful extragalactic emitters, having P> 10** watts (c/s)- m 
ster—! and p< 107! pc-3, to 
This analysis will now be extended in the light of the new observations which so 
are available. Since most of these observations were made at a frequency of de 
178 Mc/s it will be convenient to give the numerical values derived for source ch 
luminosity, volume emissivity etc. appropriate to this frequency; these quantities pr 
will differ by a factor of about two from the corresponding earlier figures. It will lu 
also be convenient, when considering extragalactic sources, to adopt the more so 
recent figure of 1tookmsec~!Mpc~ for Hubble’s constant (Sandage 1958) wl 


instead of the value 180 kmsec~! Mpc-", previously used. 

As before, a uniform but random distribution will be assumed for the region 
containing the more intense sources, witha dispersion of the luminosity of the form: 
n(P)ocexp [— (logy) P/Po)*/2] 
where n(P) represents the number of sources of luminosity P which occur in a 
sample producing flux densities within a given range. It may be rewritten to give 

the true space density p(P) of sources of different P: 

p(P) xP? exp [ — (logis P/P)?/2]. 
By analogy with an undispersed luminosity function the constant of proportion- 
ality in the first relation may be written in the form pp P,**; pp is then seen to 
represent a weighted mean density of sources. In later considerations it will be 
convenient to work in terms of pp) and P4. 

This dispersion covers a range of P similar to that used in the earlier analysis 
but its form is more probable. A smaller dispersion seems unlikely in view of the 
far larger values evident in optical stars and galaxies; the adopted dispersion 1s 
comparable with that found amongst identified ratio sources (Long and Marks 
1961). If the actual dispersion were greater than that adopted it would impose 
more stringent restrictions than those derived below. 


(6) The number-flux density relation 


At distances corresponding to small red-shift the expected number of sources 

N per unit solid angle which have a flux density greater than S is given by: 
, N=}p,P,3? S82. 

As already indicated the observations are not consistent with a uniform spatial 
density, but an observational relationship between N and S can be established at 
a large value of N where the statistical errors are small and yet not near the limit of 
the observations. The value of py Py?” derived in this way will be appropriate in 
the following analysis in which the source counts are related to the observations 
of the isotropy of the weak sources and to their contribution to the integrated 
emission. 

Using the recent observations a value of py Py? = 1-4 x 10! is found where py 
is expressed in sources pc~* and Pin watts (c/s)! ster~! at a frequency of 178 Mc/s. 
This result is shown in a plot of log P, against log py in Fig. 1. Any class of source , 
whose radio luminosity and space density are represented by a point on this line 
will satisfy the general observations of the number of sources falling in a given 
range of flux density, although not of course the non-uniformity in the observed 
relationship. 








hat 
/s) 
3 


ich 

of 
rce 
tes 
vill 
ore 


58) 


ion 


rces 


atial 
d at 
it of 
fe in 
ions 
ated 


© Po 
Ic/s. 
urce 
line 
iven 








No. 4, 1961 


(c) Isotropy 

The observations have shown that, except for low galactic latitudes and the 
area of the “‘spur’’ in the galactic emission at /" = 30°, the sources are distributed 
isotropically. As has already been demonstrated (Ryle 1958) most of the sources 
must then either be extragalactic or they must be in the galactic halo. In order 
to establish closer limits on the spatial density of the sources in the latter case, 
some new observations have been made (Scott, Ryle and Hewish 1961) which were 
designed to detect any differences in the distribution of sources in two areas of sky 
chosen to represent a large and small extent of the galactic halo. It was shown 
previously that if the smaller extent is 7, then, in the absence of dispersion in the 
luminosity, the source distribution will appear different in the two areas if the 
source density p < 3.N/r,°, where N represents the number of sources per steradian 
which can be observed. 
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Fic. 1.—Diagram relating the values of the effective radio luminosity P, and space density po 
which are compatible with the counts of sources falling in a given range of flux density. The restric- 
tions set by the observed isotropy of the sources, and the maximum contribution to the background 
emission, are also known. 
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Observations were made of a strip of sky 02°<5<08°, and two regions th 
a = 12" to 165 and «=00" to 04" were analysed by the statistical method (Scheuer es 
1957). The positions of the two areas in relation to the emitting region associated hi 
with the galactic halo (Baldwin 1955) are shown in Fig. 2. 
re 
e@ 
ar 
: 
Q 
+ ; 
Y 
Fic. 2.—The region of emission associated with the galactic halo, showing the two areas examined 
(r1=5 kpc). 
Centre of A: Mi~ 160°; Dx —55° 
Centre of B: ~355°; bU~4s°. 
The results of this analysis are shown in Fig. 3 and show no signficant differences 
in the histograms derived for the two areas. Hewish (1961) has computed the 
theoretical curves for a number of different models of source distribution and in h 
some of these a cut-off of sources below a certain flux-density was adopted to 3 
simulate the effect to be expected in the present case. Two of his theoretical tl 
curves are shown in Fig. 3 ; these correspond to similar models, but one of them has te 
a cut-off at 2 x 10 sources/steradian. The differences between the theoretical t 
curves are considerably greater than in the observed curves and it can therefore a 
be concluded that the observations indicate no anisotropy at 2 x 10* sources per ( 
steradian. 
This result shows that if the sources are mostly in the galactic halo, then with fc 
the adopted dispersion of luminosity, the effective space density py must be at $] 
least 3 x 1073 pc-, fi 
(d) The integrated emission from the sources 
Observations of the distribution of the background radiation have allowed the 7 
derivation of a model of the galactic emission, and hence the determination of the a 
average emission per unit volume in different parts of the galaxy and also in C 
extra-galactic space (Baldwin 1955). More recently, Costain (1960) has measured tl 
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the spectral distribution of some of the components and has thus been able to 
estimate the volume emissivity at different frequencies, both within the galactic 
halo and in extra-galactic space. 

Without having to know the mechanism of the emission from either of these 
regions it is clear that the contribution made by the radio sources (ZpP) cannot 
exceed, at any frequency, the volume emissivity determined in this way; this 
analysis then allows an upper limit for pp Py to be derived. 








— © Area A 
+ Ares B 
+ + (a) Model without cut- off 
(b) ° (b) Model with cut-off 
ar 
= + 
a, (a) 
z 
3 
§ 
S — 
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+ 
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l = 3 | L jl 
o 0-2 0-4 0-6 0-8 1-0 


D (equivalent fiux density x 107* w. (c/s)"' m=?) 


Fic. 3.—Histograms of recard deflections for the two areas A (anti-centre) and B (centre) 
shown in Fig. 2. The theoretical curves are also shown for two models, one of which has a cut-off 
at S=o-15 X 107" watts (c/s)-1 m=? (2 x 104 sources per steradian). 


Except for some of those at low galactic latitudes, most radio sources appear to 
have an emission which varies as Socf—8+®® over the range of frequencies 
38-960 Mc/s (Whitfield 1957; Harris and Roberts 1960). The radiation from 
the galactic halo, on the other hand, has been observed in the range 38-178 Mc/s 
to have a brightness temperature 7), ocf—(*37+°) a relation which corresponds 
to a volume emissivity proportional to f-%7+°, It is therefore clear that 
at the frequency of the present observations (178 Mc/s) only a fraction 
(38/178)©8-° of the emission from the halo can be due to the sources. 

With the adopted dispersion in source luminosity this leads to an upper limit 
for the volume emissivity of pp Py < 3:3 x 10’ watts (c/s)-!pc~*; this limit is 
shown in Fig. 1 which shows to it be incompatible with the lower limit of py derived 
from the observed isotropy of the sources. 

This analysis may be extended in the light of recent observations of the galactic 
emission which have now been made over a range of frequencies of 100: 1, and 
which can be fitted with a uniform spectral index of —0-37 (Costain 1960; 
Chapman 1960: Seeger 1960). Without the adoption of a peculiar distribution of 
the spectrum of the sources outside the range 38-178 Mc/s, it is difficult to suppose 
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that they can make a significant contribution to the observed galactic radiation; 
the combination of the radiation from two types of source having different spectral 
indices can produce a uniform spectral index only if the spectra of the two 
components show opposite departures from linearity at approximately the same 
frequency. With the uncertainties in the present observations it seems unlikely 
that the sources could contribute more than 10 per cent of the emission from the 
galactic halo at a frequency of 178 Mc/s. ‘The corresponding upper limit of 
Po Py < 6-7 x 10° watts (c/s)—? ster— pc is shown dotted in Fig. 1. 

The discrepancy between the derived limits of py Py, the lower limit of py and 
the value of py P,*!* therefore makes it difficult to suppose that more than a few 
per cent of the sources are situated within the galactic halo. 
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A similar argument may now be applied to the case in which most of the sources 
are extra-galactic since their contribution to the total volume emissivity of extra- 
galactic space must not lead to a total sky brightness greater than that observed. 
The relation between volume emissivity and sky brightness now depends on the 
cosmological model adopted (Shakeshaft 1954; Priester 1954). 


It is difficult to distinguish observationally between an extra-galactic com- 
ponent of the background emission and one due to the outermost shells of the 
galactic halo since both give rise to a distribution which is nearly isotropic when 
viewed from the position of the Sun. Owing to the differences in the spectra 
of the galactic halo and of the sources it is however, possible to determine the 
maximum contribution which sources having a given spectral distribution may 
make to the extragalactic emission. 


Costain’s observations (1960) can be used in this way to derive an upper limit 
of 25°K for the maximum sky brightness at 178 Mc/s due to sources having a 
spectral distribution Socf-**. In the case of the steady-state model and with 
the adopted dispersion in luminosity this result leads to an upper limit of 
po Py) < 150 watts (c/s)-? ster! pce-*. The intercept of this line with the line 
po Po? = 1°4 x 10" (Fig. 1) leads to a lower limit of Py) > 10% watts (c/s)! ster—. 

In most evolutionary models the total sky brightness converges more slowly 
than for the steady-state model and a lower limit must be set to py) Py; consequently 
the lower limit of the effective source luminosity P,) must be increased for such 
models. 
(e) Discussion 

The extension of the previous analysis which has been carried out using the 
new observational material has effectively eliminated the possibility that a signifi- 
cant fraction of the observed sources are situated within the Galaxy. This 
possibility is in any case hard to accept in view of the observed non-uniform 
distribution in depth; if more than a small fraction of the sources were within the 
galactic halo it would be necessary to postulate a progressive change in luminosity 
or spatial density of the sources with distance from the Sun, up to several 
kiloparsec, if both the number-flux density counts and the isotropy were to be 
explained. 

It may therefore be concluded that most of the sources are extra-galactic 
with an effective luminosity P, > 10% watts (c/s)! ster—1. 

It is interesting to compare the limiting luminosity function derived in this 
way (Fig. 4(a@)) with those derived by other arguments. 
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Measurements of the surface brightness of the sources (Allen, Palmer and 
Rowson 1960; Leslie 1961 a) allow a luminosity function to be derived if it is 
assumed that all sources are of galactic dimensions (Fig. 4(b)). An analysis 
both of sources identified with optical objects, and of the limiting optical magni- 
tudes of unidentified sources, has enabled Long and Marks (1961) to derive a 
luminosity function on the supposition that all sources have an absolute optical 
magnitude M,,= —20°6+ 1-0, a range which includes all the identified sources 
whose distances have been measured. ‘Two limiting cases of their distribution 
are shown in Fig. 4(c). In connection with their results, it is interesting to 
note that sources of intermediate radio luminosity, such as NGC 1275 and M87, 
cannot represent a large fraction of those observed in a given range of flux- 
density; some authors have suggested that the presence of such sources would 
prevent the recognition of any cosmological effect due to more powerful sources. 

Minkowski (1960), Bolton (1960) and Mills (1960) have also derived luminosity 
functions on the supposition that the sources which have been identified with 
optical objects represent a typical sample; their curves replotted to give the 
distribution in luminosity of the sources falling in a given range of flux density are 
shown in Fig. 4 (@). 
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Fic. 4.—The distribution in radio luminosity (P) of the sources falling in a given range of flux 
density. 

(a) Lower limits derived from radio observations without assumptions concerning the nature of 
the sources. 

(b) Model derived from consideration of the distribution of the surface brightness of the sources. 

(c) Model derived from limits of (m,—m,,) for identified and unidentified sources (Long and 
Marks 1961). 
(d) Models derived on basis of identified sources (Minkowski 1960; Bolton 1960; Mills 1960). 
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All the available evidence therefore indicates that most of the sources observed 
in a given range of flux density have a luminosity P of 1074-107’ watts (c/s)—'ster—!. 
Whilst arguments based on measurements of surface brightness and of optical 
luminosity are both dependent on assumptions concerning the nature of the 
unidentified sources, the results they provide agree with the lower limits which it 
has been possible to derive without such assumptions from the radio observations 
alone. 

Even if this more conservative limit is adopted, the region investigated in the 
survey of individual sources with S>2 x 10-** watts (c/s)! m~ extends to dis- 
tances of at least 10° pc. Any suggestion that the apparent non-uniformity in the 
distribution of the sources should be attributed to a local irregularity must now be 
discounted since the volume of space involved represents an appreciable fraction 
of the observable universe. 

In the next section, an investigation is made into the effects of the red-shift 
on the counts of sources having luminosities of this order. 
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3. The predicted number-flux density relationship.—Several authors have 
examined the number-flux density (NV — S) relationship expected for radio sources 
of given luminosity and have derived curves appropriate to a number of different 
cosmological models (McVittie 1957; Priester 1958; Davidson 1959). 

In most of these models evolutionary effects have been neglected and the 
density of radio sources has been supposed to vary in the same way as that of 
galaxies. The problem of accounting for the radio emission from powerful 
extra-galactic sources of the type now envisaged has however, been examined by a 
number of authors (Greenstein 1955; Burbidge and Burbidge 1957; 
Ambartsumian 1958; Burbidge 1958; Shklovsky 1960; ‘Tunmer 1960; 
Minkowski 1961). The conditions required to account for the emission involve 
a very considerable supply of energy and some of these authors have concluded that 
the kinetic energy of the gaseous content of colliding galaxies may represent the 
only adequate source—a conclusion supported by the nature of a number of the 
identified sources. 

If the radio sources have their origin in galactic collisions then their space 
density is likely to show more marked changes than that of galaxies, and in the 
limit for field galaxies would vary as the square of the mean galactic density 
(Priester 1958; Davidson 1959). 

In the present analysis two particular cosmological models have been 
considered. 


(a) The steady-state model. 
(6) The Einstein-de Sitter model. 


In the latter case evolutionary effects and the variation of the total number of 
sources may both be important. 

Curves have been plotted in Fig. 5 showing the variation with S of N/No, the 
ratio of the number of sources observed having a flux density greater than S, 
compared with the number which would be observed in a static Euclidean universe. 
In the Einstein-de Sitter model two limiting cases have been considered in which 
the density of sources is supposed to vary (i) as the density of galaxies and (ii) as 
the square of the galactic density. 
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The relations used in deriving these curves have been given by a number of 
authors (e.g. Bondi 1952; Davidson 1959). It is convenient to quote the formulae 
in terms of the red shift (1 +.) of the source and the associated distance r which is 
a function of the model adopted. pis the local density of sources and P the lumin- 
osity per steradian of the source which is assumed to obey a spectral law Pocf-*. 

For all the models considered the flux density per unit bandwidth is 

P/r?(1+2)'+* 

where P has the value appropriate to the frequency of reception, not that of the 
original emission. In the steady-state model, which has a constant number of 
sources in a given proper volume, the number of sources in the range r to (r + dr) 
is 47pr°dr/(1+2)® with r=cz/H (where H=Hubble’s constant, c=velocity of 
light). In relativistic models the number of galaxies per coordinate volume is 
constant and for the particular model chosen, the Einstein-de Sitter model, the 
corresponding relation is 4zpr’dr with r=2c/H (1—(1+2)-"?). The limiting 
case of an Einstein-de Sitter universe in which sources are produced by the 
collisions of galaxies leads to a number of sources which is proportional to the 
square of the proper density. 











l l I l 
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Fic. 5.—Curves showing the variation of the number of sources (N) having a flux density greater 
than S, compared to the number (No) on a static Euclidean model. 

(a) Steady-state (x=—o’5). 

(b) Steady-state (x= —1°o). 

(c) Einstein-de Sitter (number of sources proportional to mean galactic density). 

(d) Einstein-de Sitter (number of sources proportional to square of galactic density). 

The method of presenting the results, in which the numbers of sources are 
compared with those in a corresponding static Euclidean universe, has the 
advantage that the predictions of different models may be compared without 
a knowledge of the absolute local density of sources. 

The curves do not depend critically on the value of the spectral index x, as can 
be seen from Fig. 5, which shows two curves appropriate to the steady-state model 
with a non-dispersed luminosity P= 10% watts (c/s)-! ster-4, and x= —o°5 
and —1-o. The curves for the Einstein-de Sitter models in which the number of 
sources is assumed to be proportional (i) to the density of galaxies and (ii) to the 
square of the density, are also shown in Fig. 5; a mean spectral index x = —o-8 has 
been adopted and no evolutionary effects have been included. 
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The purpose of the present paper is not to determine a particular model which 
accounts for the observations, but to make a critical examination of the predictions 
of the steady-state model, which allows of no variations except those associated 
with the spectral index and the source luminosity function. In the next section 
the steady-state model will be compared with the observations. 


4. Comparison of the observed number-flux density relationship with that 
predicted on the steady-state model.—The number-flux density relationship at 
galactic latitudes >20° has been established directly for S>2x10-** watts 
(c/s)-! m~* by combining the ‘‘ whole-sky ’’ observations of Miss Leslie (1961 a) 
with the high resolution observations of 0-6 steradians of sky (Scott and Ryle 1961). 
By using the statistical method of record analysis on the high resolution 
observations Hewish (1961) has been able to extend the analysis to cover sources 
having 5 x 10-8 < S<2x 10-* watts (c/s)-! m~*. 

All the observations are combined in Fig. 6 which shows for S>2 x 10~%6 
watts (c/s)! m~?® the values of N/N, derived from individual source counts 
together with the limits of error due to the statistical uncertainties, possible 
effects of extended sources, source clustering etc. Three models are shown 
for S<2x10~** watts (c/s)-' m~*, which represent limiting cases of models 
compatible with the observations. 

The marked convergence shown by all three models at small S may be confirmed 
from a consideration of the contribution which the sources would make to the 
total background emission of thesky. As already indicated, the contribution to the 
extra-galactic radiation from sources having a spectral index of —o-8 cannot 
exceed about 25°K. The total emission from the sources which can be observed 
individually (S>2 x 10~** watts (c/s)-! m~*) corresponds to a sky brightness 
temperature of about 3 °K, and the variation with Sof N/N, found for S>2 x 107% 
watts (c/s)~! m~? cannot extend in the same way for S<2 x 10~*’ watts (c/s) 
m~* or the permissible integrated emission would be exceeded. The three 
models (i), (ii) and (iii) derived from the statistical analysis of the interferometric 
records have been integrated for S>10~*’ watts (c/s)-! m~ and lead to sky 
brightness temperatures of 7°5, 8-5 and 11°K, all of which are compatible with the 
observed sky brightness. 

The observed number-flux density relationship may now be compared with 
that predicted by the steady-state model, using the luminosity functions derived 
in Section 2. As can be seen from Fig. 5 the curves derived for the steady-state 
model do not depend, critically on the spectral index, and a value of — 0-8 has been 
adopted; a comparison of the intensities of a number of sources at frequencies of 
159 and 960 Mc/s (Harris and Roberts 1960) has shown that most sources have 
over this range a spectral index of 0-8 + 0-2 and errors in the predicted curve of 
N/N, are therefore likely to be small, at any rate for values of z< 5. 

Two curves have been computed and are shown in Fig. 6: the first is based on 
the limiting luminosity derived from radio evidence alone; the second on the 
assumption that all sources have the same luminosity function as that derived for 
the identified sources. The point on each curve at which the median of the 
luminosity function corresponds to z= 0-5 is shown. 

Both theoretical curves have been drawn to pass through the experimentally 
determined point at S= 30 x 10~*6 watts (c/s)"'m~*. This is the largest value of 
S for which the statistical uncertainties are small enough to provide a useful value 


of N. 
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A comparison of the predicted and observed curves shows a marked 
discrepancy, even when the smallest permissible source luminosity is adopted; 
the observed number of sources in the range 0-5 < S <2 x 10~*6 watts (c/s)-! m~* 
is 3 + 0'5 times that predicted by the steady-state model. If a luminosity function 
similar to that of the identified sources is assumed, the discrepancy is 11 + 2. 
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Fic. 6.—Comparison of the observed N|Ng relationship with that predicted by the steady-state 
model (x= —0°8). 


(a) Limiting luminosity function. 
(6) Luminosity function similar to that of the identified sources. 


The flux density at which the median value of P corresponds to z=0°5 for each of the models is 
indicated. The dotted line corresponds to the useful limit of observation in the earlier surveys 
(Bennett and Smith 1961). 


No attempt has been made here to provide an explanation for the observations 
in terms of an alternative model and it is evident that the introduction of the 
evolutionary changes which appear to be necessary will make the selection of a 
unique model difficult. These observations do, however, appear to provide 
conclusive evidence against the steady-state model. 

We should like to express our thanks to Dr A. Hewish both for useful discussions 
and for allowing us to use the results of his analysis prior to publication. One of 
us (R. W. C.) is indebted to the D.S.I.R. for a maintenance allowance. 


Mullard Radio Astronomy Observatory, 
Cavendish Laboratory, 
Cambridge: 
1961 January 28. 
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NGC 7492 








James Cuffey 
9 (Received 1961 March 10)* 
Summary 
Photometry of NGC 7492, a very faint, well resolved cluster, shows that it 
has a colour-magnitude relation similar in form to that of a typical globular 
cluster. A well defined but sparsely populated red giant branch is present, 
along with a horizontal branch, blue stars, and variables. The distance 
modulus, 17-17, and the distance, 27:2 kpc., are obtained. The colour excess 
appears to be negligible. The integrated photographic absolute magnitude is 
—4°84. The star density in the cluster is unusually low, being approximately 
0°06 stars per cubic parsec for stars brighter than + 2-5 absolute magnitude. 
nd 
NGC 7492 is a very poorly populated cluster of unusually open appearance 
well removed from the Milky Way, in southern galactic latitude — 65°. Intrinsi- 
cally, it is one of the faintest of the globular clusters. Shapley (1930) describes 
its concentration class as XII, and its angular diameter as 3'-3. Nine variable stars 
are present according to Mrs Helen Sawyer Hogg (1947), although few details 
have been published. The cluster is unusually blue judging from the colour 
excess —0-20 obtained by Stebbins and Whitford (1936). 
TABLE I 
leri. Photoelectric standards, NGC 7492 
Star P V P-V 
A 13°60 13°04 +0'56 
B 15°26 14°63 +0°63 
G 15°75 14°27 +1°48 
H 15°80 14°52 +128 
I 16°82 16°75 +0°'07 
J 16°22 15°55 +0°67 
N 16°71 15°88 +0°83 
Oo 17°79 16°87 +0°'92 
P 17°93 18-26 — 0°33 
Q 16°43 15°73 +0°70 
R 16°43 15°63 +080 
S 16°11 15°54 +0°57 
. 15°81 15°16 +0°65 
U 16°77 15°78 +0'99 
g 17°08 16°32 +0°'76 


The photometric measures of cluster stars reported here are based on 
photoelectric standard stars observed in 1951 with the 82-inch reflector at the 
McDonald Observatory (Table I), and on iris photometer measures of 








* Received in original form 1960 December 28. 


photographs taken with the 1oo-inch and the 200-inch reflectors of the 
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Mt Wilson and Palomar Observatories (Table II). A dry ice refrigerated 
1P21 photometer was employed in the photoelectric work, and the measures 
were later reduced to the PV system by measures with the 36-inch reflector 
of the Goethe Link Observatory to relate them to the PV (Eggen) and the 
BV (Johnson and Morgan) systems. ‘These reductions appear in another 
publication (Cuffey and McCuskey 1956), and are very linear in nature. The 
values in Table I are reduced to the PV system (Eggen 1955). 


TaB_e II 
NGC 7492, photographic observational material 
Exposure 
Object Plate no. (min.) Colour Instrument Aperture 
NGC 7492 Hs © 30 B 100 58 
PH 32 O 10 B 200 200 
H 109 C 60 B 100 58 
H 119 C 60 B 100 58 
444 C 30 B 100 58 
NGC 7492 H 320 15 ¥ 100 58 
PH 33 O 20 = 200 200 
H 110 C 60 ¥ 100 58 
H 120 C 60 ¥ 100 58 
H 133 C 30 ¥ 100 58 


The photographic observations (Table I1) made at the Mt Wilson and 
Palomar Observatories employed 103aO plates with a WG2 filter, and 103aD 
plates with a GG14 filter to cover the photographic and the yellow regions, 
respectively. ‘Iwo of the plates, taken with the 200-inch reflector, and two of 
those taken with the 100-inch telescope were made by Dr Donald Osterbrock, 
who had begun a study of the cluster, and who very kindly presented them to me. 
The rest of the plates were taken by the author in 1955, as guest investigator at 
Mt Wilson Observatory. In addition to the plates listed, several photographs 
of Selected Areas 57 and 68 were made, and a number of similar plates were 
borrowed from Dr A. R. Sandage. 

The photoelectric observations with the 82-inch reflector gave the zero point 
and scales of the P and V magnitudes in Table I in the range 12:0 to 18-4. 
However, in order to provide more numerous standard stars for the photographic 
reduction curves, supplementary stars were added, and their magnitudes 
determined by photographic interpolation and extrapolation. In order to 
perform this work, it was essential to determine accurately the form of the iris 
reading—magnitude relation, and the photographs of Selected Areas 57 and 68 
were measured for this purpose. Photoelectric magnitudes and colours in the 
Selected Areas (Stebbins, Whitford and Johnson, 1950) were used, and also 
values by W. A. Baum, who kindly permitted me to use his photoelectric values 
in Selected Areas 57 and 68 in advance of publication. The Selected Area plates 
established the form of the iris photometer reduction curves, and showed them 
to be very smooth and free from bumps or wiggles. Therefore they could be used 
for accurate interpolation between the photoelectric standards. ‘They also 
established the form of the reduction curves for the fainter stars so that 
extrapolation in the range 18-4 to 19-5 could be accomplished, although, it must 
be admitted, with considerably less confidence than the interpolation. 
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Table III summarizes the magnitudes and colours of all the standard stars, 
both those obtained photographically by interpolation and extrapolation of the 
photoelectric data, and also those re-derived for the stars in Table I. The zero 
point was derived from photoelectric comparisons with the North Polar 
Sequence, made with the 82-inch, and using the revised values of magnitudes 
which represent the P and V system, rather than with the older International 
Polar Sequence values. The scale is that of the photoelectric measures with the 


Tasxe III 
NGC 7492, standard stars 


Star to V P-V 

A 13°77 13°13 0°64 

B 15°14 14°55 0°59 

Cc 16°90 16°71 o'19 

D 16°82 16°37 0°45 

E 16°99 16°30 0°69 

F 17°20 16°80 0°40 

G 15°78 14°18 1°60 

H 15°79 14°49 1°30 

I sea wae co D 
J 16°37 15°57 0°80 

K 16°31 15°49 082 

L 16°58 15°82 0°76 

M 16°86 16°16 0°70 

N 16°85 16°07 0°78 

O 16°81 16°07 0°74 

% 17°18 17°00 o'18 V 
Q 16°60 15°77 0°83 

R 16°25 15°31 0°94 

i) wii ee a D 
é 16°09: 15°24: 0°85: D 
U 16°59 15°81 0:78 

V 16°69 16°15 0°54 

W 17°05 16°72 0°33 

X 17°45 16-98 0°47 

¥ 17°41 17°56 —O'ls 

Z 17°17 16°79 0°38 

a 17°35 17°50 O85 

b 17°09 16°43 066 

c 17°43 16°74 0°69 

d 17°56 16°99 0°57 

e 16°82 16°24 0°58 

f 17°06 17°23 —0O'17 

g 16°89 16°19 0'70 

h 17°22 16°65 0°57 

i 17°27 17°44 —0'17 


D=double star, V = variable. 


82-inch reflector, and represents the photometric system of the present paper. 
The values could undoubtedly be improved by further measurements with a 
sensitive cell and a larger telescope, especially for the fainter stars, but they 
may be considered satisfactory for the present preliminary study of the cluster. 

Magnitudes and colours for the resolved stars in NGC 7492 are given in 
Table IV, and the identifications are given by the chart.in Fig, 1. Probable 
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S1R2 
17°25 16°85 
17°05 17°23 
20°21 19°34 
19°92 18°94 
19°55 19°16 
19'29 18°59 
20°69 19°44 
17°50 17°74 
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20°09 18-91 
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17°51 17°66 
17°38 17°55 
17°53 16°96 
17°09 17°22 
16°56 15°76 
19°04 18°20 
19°57 18°73 
17°47 16°85 


S1R3 
19°84 18-92 
18-28 17°25 
19:26 17°78 
20°44 19°37 
18-o1 17°45 
19°63 18°43 


S2R1 
18°26 17°75 
19°83 17°27 
16:90 16°40 
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+0°40 
—o'18 
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+096 
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—O'ls 
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Taste IV 
NGC 7492, catalogue of stars 


V 


18:21 
18-62 
17°37 
16°89 
17°80 
17°34 
18°47 
19°07 
16°90 
17°62 
17°71 
16°39 
16°15 
17°61 
18°52 
18-96 
19°13 


S2R2 


19‘10 
19°49 
19°14 
20°19 
18°81 
17°38 
Ig‘II 
17°28 
18-86 
19°22 
19°45 


18-32 
18°82 
18°53 
19°29 
18-10 
17°56 
18-38 
17°48 
18-17 
18°73 
18°63 


S2R3 


18°80 
19°42 
15‘10 
18-38 
20°01 


18-16 
18°61 
14°53 
17°35 
19°22 


S$3R1 


17°07 
18-99 
18-62 
18-92 
18°40 
16°95 
17°38 
16°70 
17°11 
18-22 
18°95 
18°77 
19°00 
17°44 
17°46 





17°27 
18-32 
18-00 
18°25 
17°82 
16-92 
16°77 
16:00 
17°25 
17°79 
18°33 
18-30 
18-27 
17°56 
17°58 





P-V 
+0°74 
+0°86 
—0°17 
+0°56 
+0°48 
—o°l4 
+0'50 
+0°73 
+0°53 
—0'09 
+0°53 
+0°66 
+0°57 
—O'l4 
+0°54 
+ 1°04 
+0°'90 


+0°78 
+0°67 
+0o0°61 
+090 
+o°7I 
—o'18 
+0°73 
—0°20 
+0°69 
+0°49 
+082 


+0°64 
+o'81 
+0°57 
+ 1°03 
+0°79 


—0'20 
+0°67 
+0°62 
+0°67 
+0°58 
+0'03 
+0o°61 
+0°70 
—o'l4 
+0°43 
+0°62 
+0°47 
+0°73 
—Oo'12 
—o'l2 
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17°53 
19°25 
17°46 
17°94 


V 
17°62 


18°54 
16°81 


17°31 


S3R2 


20°52 
19‘I2 
19‘OI 
19°77 
20°18 
18°57 
17°25 
18-09 
17°34 
16°68 
16°80 


19°45 
18°37 
18°31 
18-66 
18-92 
17°89 
16°60 
17°53 
17°45 
15'92 
15°99 


S3R3 


19°25 
1g‘I2 
19°43 
20°16 
20°30 
13°75 
18°83 
18°73 


18°34 
18-38 
18°63 
18-96 
19°00 
12°86 
18°07 
17°88 


S4R1 


18°12 
19°36 
18°35 
18-18 
17°14 
16°82 
18-66 
19°28 
16°94 
17°11 
17°17 
17°41 
16-96 
18°03 


17°48 
17°68 
17°77 
17°63 
17°18 
16°21 
18-09 
18-38 
17°03 
17°28 
17°28 
16°77 
16°23 
17°41 


S4R2 


17°58 
18-89 
18-86 
18-65 
18-32 
18°51 
20°34 
17°23 
19°14 


17°70 
18-19 
18°14 
17°99 
17°67 
17°99 
19°05 
17°32 
18°43 
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P-V 
—0'09 
+0°71 
+0°65 
+0°63 


+1°07 
+0°75 
+0°70 
+ III 
+1°26 
+0°68 
+0°65 
+0°56 
—o'll 
+0°76 
+081 


+091 
+0°74 
+080 
+ 1°20 
+ 1°30 
+0°89 
+0°76 
+0°85 


+0°64 
+1°68 
+058 
+0°55 
— 0°04 
+0°61 
+0°57 
+0°90 
—0-09 
—0O'17 
-o'll 
+0°64 
+0°73 
+0°62 


—O'12 
+0°70 
+0°72 
+066 
+0°65 
+0°5§2 
+ 1°29 
—-0°09 
+0°71 
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TaBLe IV (cont.) 
No. | a V P-V No. P V P-V No. P V P-V 
S4R2 16 18:01 17°50 +0°51 3 16°30 14°83 +1°47 
IO 17°23 17°31 —0°'08 17 18:90 18°14 +0°76 4 20°39 19°00 +1°39 
It 18:23 17°68 +0°55 18 19°67 18°81 +0°86 5 20°39 18°94 +1°45 
12 18:52 17°84 +0°68 6 17°70 17°10 +0°60 
13 18°68 18-08 +060 S4R3 7 19°65 18:92 +0°73 
14 18-81 18:24 +0°57 I 20°22 19°00 +1°22 8 20°47 19°21 +1°26 
1§ 17°24 17°37 —0°13 2 18:15 17°60 +0°55 9 20°02 18°87 +115 


errors, which reflect only the internal consistency of the data, are summarized 
as follows: brighter than V=18-5, +0-04 in V, and +0:06 in P—V; fainter 
than V=18-5, +0:07in V and +o:10in P—V. There are several factors which 
rendered the present photometry of NGC 7492 unusually difficult from the 
standpoint of obtaining the highest accuracy. The cluster’s declination is 
16° south, and the Mt Wilson plates, taken near the meridian, encountered 
somewhat higher fog than usual from the brightly illuminated sky over Pasadena. 
Improved photometry could result from observations made at a southern 
hemisphere observatory. The photoelectric observations, made at McDonald 
with the 82-inch reflector, however, were made under good observing conditions, 
and the magnitudes of Table IV are based on photoelectric standards to magnitude 
V = 18-4, so that they are probably unaffected by plate fog uncertainties to this 
limit. Fainter than V=18-o0, the present results should be regarded as pro- 
visional. 

The colour—magnitude relation, Fig. 2, shows that in spite of its poverty 
in stars and its open structure, NGC 7492 is a globular cluster. The red giant 
branch is well defined, and stars as red as P— V = 1-6 lie on it. It is, however, 
poorly populated towards the red end, and this fact, along with the well populated 
horizontal branch, probably explains the negative colour excess obtained by 
Stebbins and Whitford (1936) who derived the colour excess by comparing 
NGC 7492 with other globular clusters. The horizontal branch reaches 
P—V = -—0o-24, and this, combined with the cluster’s high south galactic latitude 
(—65°), makes it probable that selective absorption is negligible. Fainter than 
V =18-0, the observations have been plotted as open circles to indicate their 
tentative nature. 

Nine variable stars are present, according to Mrs Helen Sawyer Hogg (1947). 
Much of the material, however, is unpublished, and insufficient data are now 
available on the variables to make them useful for distance determination in 
NGC 7492. Shapley (1920) listed one variable and four suspects. The variable, 
V1, is standard star P and it shows large residuals in the measures, thus confirming 
its variability. The suspects are, in the designation system of Table IV, S1R2 19, 
standard star f, and S4R21; Shapley suspect number 5 lies outside the area 
measured. None of the suspects investigated showed unusually large residuals 
in the present series of measures. Furthermore, only V1 lies at the expected 
location in the colour-magnitude relation; suspect V2 lies too far to the red 
to be on the horizontal branch, while V3 and V4 are too blue. Stars which 
should be further investigated for variability because of their location on the 
horizontal branch of the colour-magnitude relation are standard stars C, F, W, 
and Z, and stars S1R17, S1R21, S3R16, S4R15, and S4R19. The present 
series of plates is insufficient to give further information on the variables. 
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Fic. 1.—NGC 7492 star chart. Standard stars are circled. 
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The main sequence lies below the observational limits of Fig. 2. Fainter 
than V = 18-0, the measures become uncertain because of interfering plate fog 
from the bright sky over Pasadena, and it is not possible to be sure that the sharp 
trend toward red stars at V=1g:0 is real and not caused by systematic errors 
in the colours. Until a photoelectric check of the sequences at V= 190 can be 
made, little attention should be paid to trends indicated in Fig. 2 fainter than 
V = 18-0. 

















Fic. 3.—Comparison of colour-magnitude relations in NGC 7492 and Messier 3. 


The distance for NGC 7492 was obtained by comparing its colour-magnitude 
relation with that of Messier 3 (Johnson and Sandage 1956). A zero-point 
shift of o-17 magnitudes was used to allow for the difference between the P—V 
and B—V photometric systems. Since evidence for selective absorption is 
slight in both cases, no allowance was made for it. For Messier 3, the modulus 
15°67 was used (Sandage 1953). From the blue stars the modulus for NGC 7492, 
17°24, was obtained; the horizontal branch gives 16-91, and the red giant branch 
gives 17°37. We adopt the mean, or 17-17, as the best value, and the distance is 
27'2kpe. The integrated absolute photographic magnitude is then — 4:84, if 
the apparent magnitude 12-33 (Christie 1940) is used. 

The cluster’s diameter, estimated from the Palomar Atlas prints, is 3'-3 if 
the nucleus only is considered, and 4’-5 if the outlying stars are included. The 
nucleus is therefore 26-0 pc in diameter, and the whole cluster, 35-6 pc. The 
cluster’s total population, brighter than V=19°5, was counted as 575 stars. 
Its mean density, for stars brighter than + 2:5, is 0-02 stars per cubic parsec; 
its central density, obtained from an analysis of star counts by Plummer’s method, 
is 0°06 stars per cubic parsec. 

The star counts showed that 8 per cent of the points in Fig. 2 can be attributed 
to field stars, so that the scatter in the colours of the fainter stars can be more 
probably blamed upon photometric uncertainties than interfering field stars. 

NGC 7492 is located one full Milky Way diameter from the Milky Way 
plane, in fact, somewhat farther from it than the Magellanic Clouds. Two other 
clusters, NGC 5053, and 5466 resemble NGC 7492 in structure and in location 
well removed from the Milky Way plane. If they represent objects populating 
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intergalactic space, and are unrelated to the galactic system, it is tempting to 
speculate that they owe their common poverty in stars to the scarcity of star 
building material in intergalactic space at the time of the cluster’s origin. 
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